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$ 

incompressible  materials 

by 

Rohan  C.  Abeyaratne 
California  Institute  of  Technology 

Summary 

This  investigation  is  concerned  with  the  possibility  of  the  change  of 
type  of  the  differential  equations  governing  finite  plane  elastostatics  for 
incompressible  elastic  materials,  and  the  related  issue  of  the  existence  of 
equilibrium  fields  with  discontinuous  deformation  gradients.  Explicit 
necessary  and  sufficient  conditions  on  the  deformation  invariants  and  the 
material  for  the  ellipticity  of  the  plane  displacement  equations  of  equilibrium 
are  established.  The  issue  of  the  existence,  locally,  of  "elastostatic  shocks"— 
elastostatic  fields  with  continuous  displacements  and  discontinuous  deforma- 
tion gradients  — is  then  investigated.  It  is  shown  that  an  elastostatic  shock 
exists  only  if  the  governing  field  equations  suffer  a loss  of  ellipticity  at  some 
deformation.  Conversely,  if  the  governing  field  equations  have  lost  ellip- 
ticity at  a given  deformation  at  some  point,  an  elastostatic  shock  can  exist, 
locally,  at  that  point.  The  results  obtained  are  valid  for  an  arbitrary  homo- 
geneous, isotropic,  incompressible,  elastic  material. 


< The  results  communicated  in  this  paper  were  obtained  in  the  course  of  an 

investigation  supported  by  Contract  N00014-75-C-0196  with  the  Office  of 
Naval  Research  in  Washington  D.  C. 
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1 . 1 Introduction 

In  two  recent  papers  [l],  [2],  Knowles  and  Sternberg  looked  in- 
to the  question  of  the  possible  loss  of  ellipticity  of  the  displacement  equa- 
tions of  equilibrium  of  nonlinear  elastostatics  for  compressible  materials. 
In  [l],  three  dimensional  homogeneous  deformations  of  a particular  iso- 
tropic compressible  elastic  material  were  considered,  and  necessary 
and  sufficient  restrictions  on  the  principal  stretches  for  ellipticity 
to  prevail  were  deduced.  It  was  shown  that  for  this  material,  a loss  of 
ellipticity  occurred  at  sufficiently  severe  local  deformations.  In  [2] 
they  established  similar  explicit  necessary  and  sufficient  conditions  for 
an  arbitrary  homogeneous,  isotropic,  compressible  elastic  solid  sub- 
i jected  to  plane  deformations. 

These  papers  were  motivated  by  some  asymptotic  studies  of 
crack  problems  they  had  considered  previously,  in  which  certain  diffi- 
culties encountered  suggested  that  the  problem  may  not  admit  a classi- 
cally smooth  solution. 

In  a subsequent  paper  [3]  Knowles  and  Sternberg  investigated 
the  implications  of  a loss  of  ellipticity.  The  question  of  the  existence 
of  "elastostatic  (or  equilibrium)  shocks"  — solutions  which  possess  fi- 
nite jump  discontinuities  of  the  displacement  gradient  across  certain 
surfaces  while  maintaining  continuous  displacements  — was  studied 
• within  the  context  of  plane  deformations  of  compressible  elastic  solids. 

It  was  established  in  [3]  that  a necessary  condition  for  the  existence  of 
a piecewise  homogeneous  elastostatic  shock  was  that  the  material  lose 
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strong  ellipticity  at  some  homogeneous  deformation.  The  question  of 
whether  in  fact  a loss  of  ordinary  ellipticity  was  necessary  was  left 
unanswered.  In  the  particular  case  of  weak  elastostatic  shocks  it  was 
shown  that  ordinary  ellipticity  must  necessarily  be  lost  at  the  pre- 
assigned deformation  on  one  side  of  the  shock. 

Rice  [5]  had  previously  examined  the  phenomenon  of  "localiza- 
tion of  deformation"  for  plastic  materials.  Localization  is  the  bifurca- 
tion of  an  initially  smooth  state  of  deformation  into  one  involving  a zone 
of  highly  localized  shearing.  Localized  deformations  as  described  in 
[5]  appear  to  have  certain  qualitative  features  in  common  with  elasto- 
static shocks  as  described  in  [3].  In  fact,  within  his  setting.  Rice 
shows  that  the  onset  of  localization  is  first  possible,  in  a program  of 
deformation,  when  the  displacement  equations  of  equilibrium  lose  el- 
lipticity. 

In  the  present  study  we  treat  the  corresponding  issues  for  an 
arbitrary  homogeneous  incompressible  elastic  solid  subjected  to  plane 
deformations.  Some  of  the  results  established  are  appropriate  only 
for  isotropic  materials.  Explicit  necessary  and  sufficient  restrictions 
on  the  deformation  invariants  and  the  material  are  deduced  which  ensure 
ellipticity  of  the  plane  displacement  equations  of  equilibrium.  In  the  con- 
text of  isotropic  materials  it  is  established  that  a loss  of  ordinary  ellip- 
ticity at  some  homogeneous  deformation  is  a necessary  condition  for  the 
existence  of  a piecewise  homogeneous  elastostatic  shock.  It  is  further 
shown  that  a strict  loss  of  ordinary  ellipticity  at  a given  homogeneous 
deformation  is  sufficient,  but  not  necessary,  to  ensure  the  existence  of 
a piecewise  homogeneous  elastostatic  shock  which  has  associated  with 
it  this  pre  as  signed  deformation  on  one  side. 


In  Section  2 we  cite  some  relevant  results  from  the  theory  of  fi- 
nite elastostatics  for  incompressible  elastic  solids  which  we  then  spe- 
cialize to  plane  deformations.  The  notion  of  the  "local  amount  of  shear" 
associated  with  any  plane  volume  preserving  deformation  is  then  de- 
scribed. In  Section  3 the  conventional  notion  of  ellipticity  is  adapted 
to  the  displacement  equations  of  equilibrium  and  necessary  and  suffi- 
cient conditions  for  ellipticity  are  then  deduced.  In  the  isotropic  case 
these  conditions  are  put  into  explicit  form  and  a simple  interpretation 
is  given  in  terms  of  what  we  call  the  "local  amount  of  shear".  A loss 
of  ellipticity  is  found  to  depend  on  a loss  of  invertibility  of  the  shear 
stress -amount  of  shear  relation  in  simple  shear. 

The  notion  of  piecewise  homogeneous  elastostatic  shocks  devel- 
oped in  [3]  for  the  compressible  case  is  extended  to  the  incompressible 
case  in  Section  4.  In  Section  5 we  then  consider  weak  elastostatic  shocks 
in  homogeneous,  incompressible,  anisotropic  elastic  solids  and  show 
that  a loss  of  ellipticity  at  the  pre-assigned  deformation  on  one  side  of 
the  shock  is  necessary  for  its  existence.  The  jumps  across  the  shock- 
line of  various  physically  significant  field  quantities  are  also  deduced. 

In  Section  6 we  return  to  equilibrium  shocks  of  finite  strength  in  homo- 
geneous, incompressible,  isotropic  elastic  solids.  We  show  that  a strict 
failure  of  ordinary  ellipticity  at  a given  deformation  is  sufficient  to  ensure  the 
existence  of  a piecewise  homogeneous  elastostatic  shock  which  has  associated 
with  it  this  deformation  on  one  side.  Moreover  we  show  that  a failure  of  ordinary 
ellipticity  at  some  homogeneous  deformation  is  necessary  for  the  existence 
of  a shock  of  the  type  under  consideration. 

In  Section  7 we  discuss  the  dissipativity  inequality  first  proposed 
by  Knowles  and  Sternberg  in  { 3 ] and  later  extended  by  Knowles  [4]  to 
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three- dimensional  deformations  of  both  compressible  and  incompressible 
materials  and  explore  some  of  its  consequences.  In  particular  its  im- 
plications in  the  case  of  weak  elastostatic  shocks  in  anisotropic  materials 
is  examined. 

Finally  in  Section  8 we  illustrate  some  of  the  preceeding  results 
by  means  of  an  example  involving  a particular  hypothetical  constitutive 
law. 


2.  1 Preliminaries  on  Finite  Plane  Elastostatics 

Let  R be  the  three-dimensional  open  region  occupied  by  the 
interior  of  a body  in  its  undeformed  configuration.  A deformation  of  the 
body  is  described  by  a sufficiently  smooth  and  invertible  transformation 

y = y(x)  = x + u(x)  on  R (2.1) 

which  maps  R onto  a domain  RA  . Here  y is  the  position  vector  after 
deformation  of  the  particle  which,  in  the  undeformed  configuration  was 
located  at  x . We  will  assume  for  the  moment  that  the  displacement 
vector  field  u(x)  is  twice  continuously  differentiable  on  R . 

The  deformation  gradient  tensor  F is  defined  by 

F = on  R , (2.2) 

and  since  the  material  is  presumed  to  be  incompressible, 

det  F = 1 on  R , (2.3) 

where  det  F is  the  Jacobian  of  the  mapping  (2.  1).  The  right  and  left 
Cauchy-Green  tensors  C and  G are  defined  respectively  by 

C = FTF  , G = FFT  . (2.4) 

<"W  »v 

Let  jr  be  the  Cauchy  stress  tensor  field  accompanying  the  de- 
formation at  hand.  Assuming  that  jt  is  continuously  differentiable  on 
R^  , the  equilibrium  equations  are 


div  t = 0 , T = T on  R*  , (2.5) 

rs*  /v*  T ' ' 

where  body  forces  are  presumed  to  be  absent.  The  nominal  (Piola) 
stress  tensor  corresponding  to  r is  given  by 


T - 1 

0 = t(FX)  , 


(2.6) 


where  use  has  been  made  of  (2.3).  Equations  (2.2),  (2.3),  (2.5)  and 
(2.6)  lead  to  the  equilibrium  equations 


div  a = 0 , oFT  = FoT  on  ft  . 


(2.7) 


We  now  turn  to  the  constitutive  relations  and  suppose  that  the 
body  is  homogeneous  and  elastic  and  possesses  an  elastic  potential 

A 

W = W(F)  . W represents  the  strain  energy  density  per  unit  undeformed 
volume.  The  nominal  stresses  are  now  given  by 


o = WF(F).p(FT)-1  , 


(2.8)' 


where  p(x)  is  a scalar  field  arising  because  of  the  incompressibility 
constraint.  We  assume  for  the  moment  that  p(x)  is  continuously  dif- 
ferentiable on  ft  . Alternatively,  from  (2.6),  (2.8)  we  have 


t = WT-,(F)F1  - pi 


From  (2.  l)-(2, 3),  (2.  7)  and  (2.  8)  we  are  led  to 


(2.9) 


*See  Truesdell  and  Noll  [6],  page  304. 
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c...  .(F)u.  ..-p  .F..*  = 0 

~ k,  <3  ji 


on  ft 


(2. 10)1 


where 


82W(£) 

9F..9F,  , ’ 
ij  kl 


(2. 11) 


2 2 2 

Let  X j(x)  , X2(x)  and  X3(x)  * where  XjX)  , be  the  eigenvalues 
of  the  symmetric  positive  definite  tensor  field  G (or  C).  The  princi- 
pal scalar  invariants  of  G are 


I1=  tr  G =X  2 +X  2 +X  2 . I2=-|[  (tr  G)2-  (tr  G2)]  =X  ^X  2 +X  ^X2  +X  ^X  2 , 

I3=detG  = X2X^X2  , (2.12) 

where  tr  denotes  the  trace.  From  (2.3),  (2.4)  and  (2.  12)  it  follows 
that 


on  ft  . 


(2.  13) 


In  the  special  case  when  the  material  is  incompressible  and  isotropic, 
W depends  on  F only  through  the  invariants  Ij  and  I2  , whence 
we  have 


All  tensor  and  vector  components  are  taken  with  respect  to  a fixed 
rectangular  cartesian  frame.  A comma  followed  by  a subscript  indi- 
cates differentiation  with  respect  to  the  corresponding  x-coordinate. 
Latin  subscripts  take  the  values  1,2,3  while  Greek  subscripts  take  the 
values  1,2.  Repeated  subscripts  are  summed  over  the  appropriate 
range. 
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W = W(Ij  , i2) 


(2.14) 


Suppose  now  that  the  domain  R occupied  by  the  undeformed 
body  is  a right  cylinder  with  generators  parallel  to  the  x^-axis.  Let 
II  be  the  open  region  of  the  (x^ , x^J-plane  occupied  by  the  interior  of 
the  middle  cross-section  of  this  cylinder.  Suppose  further  that  the 
deformation  (2.1)  is  a plane  deformation  so  that 


ya=  VUo(xl'x2)  ' y3  = X3  on  * 


(2.  15) 


II  is  then  mapped  onto  a domain  TI^  of  the  same  plane,  which  would  be 
the  middle  cross-section  of  the  cylindrical  region  R^  . From  here  on 
we  shall  be  exclusively  concerned  with  plane  deformations  unless  spe- 
cifically stated  otherwise.  It  follows  from  (2.2)  and  (2.  15)  that 


Fc*|3  ya,  (3  ' Fa3"F3oT0  ' F33"  1 


The  nominal  stresses  are  now  given  by 


8W(F)  _j  9W-'F) 

°ap=  " PFpo  ' °33=  dF^~  - P 


If  we  assume  that  the  elastic  potential  W is  such  that 


(2.16) 


(2.17) 


■ • 


9W(F)  eW(F) 

9F  . = 9F,  ~ = 0 (2.18) 

ct3  3a 


for  every  F such  that  (2.  16)  holds,  then  we  further  have 


= 0 . 


(2.  19) 


The  assumption  (2.  18)  holds  true  identically  for  isotropic  materials  in 
particular. 

One  sees  readily  from  (2.7),  (2.15)-(2.  19’  that  for  equilibrium 
in  the  -direction  it  is  necessary  and  sufficient  that  the  scalar  field 
p(x)  be  independent  of  . Thus 

p = p(Xj,x2)  on  II  . (2.20) 


In  the  present  circumstances  (2.  10)  specializes  to 

CaPY6(~)Uv,p6‘p,  pFpo  = 0 on  11  * (2.21) 

Equation  (2.21),  together  with  the  incompressibility  condition  (2.3), 
constitute  the  governing  system  of  equations  for  the  plane  strain 
problem  and  we  shall  refer  to  them  as  the  displacement  equations  of 
equilibrium  in  plane  strain.  They  are  three  scalar  equations  involving 
the  three  functions  u^fx^x^)  p(Xj,x2)  * 

One  sees  readily  from  (2.4)  and  (2.  15)  that  in  any  plane  defor- 
mation, unity  is  an  eigenvalue  of  the  left  and  right  Cauchy-Green  ten- 
sors, whence  we  have 
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Equations  (2.  12)  and  (2.  13)  now  specialize  to 


I1  = X1+X2  + 1 » I2=X1  + X2  + X 1X2  » r3  = XlX2  » 


and 


x2x2-i 
K 1A  2 " 1 ' 


whence 


I1-I2=X  1 + 2 + 1 
X 1 


If  we  now  define  I by 


(2.23) 


(2.24) 


(2.25) 


I=FopFaP  ' 


we  have,  because  of  (2.4),  (2.12),  (2.16)  and  (2.25)  that 


(2.26) 


I=I1-1  = I2-1  = 2 + 


*2  . 


(2.27) 


In  the  special  case  when  the  material  is  isotropic,  we  have 
from  (2.  14)  and  (2.27)  that,  in  plane  deformations, 

$ 

W = W(I+  1 , 1+  1) 


(2.29) 


so  that  if  we  define  the  Plane  Strain  Elastic  Potential  W(I)  by 
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W(I)  = W(I+  1,1+1),  I*  2 , 


(2.30) 


rN 

we  have  in  the  present  context  that  W(F)  = W(I)  where  I = . 


It  follows  from  this  that 


8W(F) 

= V1'  • 


(2.31) 


2 - 

a w(F) 

coPv6(£)=  8F;-p5V-6  = 2^vVw  (1)  + 4Vv*w  (I)  • 


(2.32) 


From  (2.4),  (2.9)  and  (2,31)  we  conclude  that 


>=2W/<I)Gap-P6ap  * 


(2.33) 


It  is  apparent  that  the  plane  strain  elastic  potential  W(I)  fully  determines 
the  in-plane  stress  components.  This  is  not  true,  however,  of  the  com- 
ponent . 

Finally  we  recall  that  in  this  case  the  in-plane  Baker-Ericksen 
inequality  requires  that 


<T  1“  T2)(x  rX2)>0  lf  X 1^X2 


(2. 34)" 


for  all  pure  homogeneous  (plane)  deformations  of  the  form 


y =X  x (no  sum);  X .X  1 , X >0  , y.  = x..  , 

QrQfQf  A *■  J j 


(2.35) 


See  Truesdell  and  Noll  [6],  page  158. 
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where  are  the  principal  in-plane  Cauchy  stresses.  From  (2.4), 

(2.  16),  (2.  26),  (2. 33)  and  (2. 35)  we  have 

Ta=2W'(I)x£-p  , I = X^+X^  , (2.36) 

whence  (2.  34)  may  be  equivalently  written  as 

W'(X  j + X^HX  j-X2)2(X  j + X2)>0  , X 1,X2>0  , X jX2=  1,  X l1t\2  (2.37) 

or 

W#(X  j +X  2>>0  , XrX2>0  , X jX  2=  1 , X j f<X  2 , (2.38) 

which  in  turn  is  equivalent  to 

W'(I)>0  for  I>2  . (2.39) 

The  infinitesimal  shear  modulus  is  easily  shown  to  be  & =2W'(2)  ; if 
we  assume  that  &>0  , we  may  replace  (2.39)  by 

W'(I)>0  for  I*  2 . (2.40) 

Requiring  that  (2.40)  hold  for  the  material  at  hand  is  equivalent  to  re- 
quiring that  the  mat-srial  have  a positive  (finite)  shear  modulus.  Con- 
versely, (2.40)  implies  (2.34),  though  it  does  not  imply  the  full  (three- 
dimensional)  Bake r- Ericksen  Inequalities. 

2.  2 Local  Amount  of  Shear 

We  now  establish  that  any  plane  volume  preserving  deformation 
can  be  decomposed  locally  into  the  product  of  a simple  shear  in  a suit- 
able direction  followed  or  preceded  by  a suitable  rotation. 
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To  this  end,  let  F be  a two-dimensional  tensor  such  that 
det  F = 1.  Define 

k =\/l  - 2 , IrF^pF^p  . (2.41)1 

Then  we  will  show  that  there  exist  proper  orthogonal  tensors  Qj  , , 

non-singular  tensors  Kj  , with  unit  determinant  (all  two-dimen- 
sional) and  rectangular  cartesian  frames  Xj  , X2  such  that 

~ = ~1  ~1  = ~2^2  ’ (2.42) 


(2.43)2 


Conversely,  if  (2.42)  holds  for  some  proper  orthogonal  tensors  Qj, 

and  tensors  Kj,  K 2 with  unit  determinant  such  that  (2.43)  is  true  in 

some  rectangular  cartesian  frames  Xj,  X2  , then  we  will  show  that 

k is  necessarily  given  by  k = ±Jl  - 2 . 

In  order  to  prove  the  first  part  of  the  result,  let  X be  a princi- 

T 

pal  frame  for  the  symmetric  positive  definite  tensor  FF  . Then 


(FFT)X  = 


-2 


\ >0 


(2.44) 


*Since  det  F = 1 , we  have  that  necessarily  I*  2. 

2 X1 

K,  1 is  the  matrix  of  components  of  the  tensor  K.  in  the  frame  X.  . 

^ 1 ~ 1 *w| 
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where  we  have  made  use  of  the  fact  that  det  F = 1.  Clearly  we  may 
assume  X £ 1 with  no  loss  of  generality.  Consider  the  rectangular 
cartesian  coordinate  frame  X^  obtained  by  rotating  the  frame  X 
counterclockwise  through  an  angle  8 determined  by 


sin  8 = 


cos  8 = 


By  the  change  of  frame  formula  for  tensors. 


(2.45) 


T T X t 

(FF1)  = R(FF  1 ) R1  , R = 


cos  8 sin  6 


-sin  8 cos  6 


T ^2 

we  compute  (FF  ) to  find 


x2+x'2- 1 X -x_1 
x-x'1  1 

Let  be  the  tensor  with  unit  determinant  defined  by 

\\  x-x'1! 


0 1 


T *2 
(FFA)  * 


(2.46) 


(2.47) 


Then  (2.46)  and  (2.47)  imply  that 


fft=k?kJ  . 


(2.48) 


Define  the  tensor  Q0  by 
~c 


Q 


2 


(2.49) 


(2.48)  and  (2.49)  now  lead  to 


T 

Since  K,  is  non- singular  it  thus  follows  that  Q0Q_  = 1 whence  Q_ 

is  orthogonal.  But,  from  (2.  49)  it  follows  that  det  Q2  = + 1 since 

detK2=detF  = 1 , so  that  in  fact  is  proper  orthogonal. 

Finally,  since  we  are  assuming  X ^ 1 , it  follows  from  (2.46) 

that  X - X * = VF^pF^p-  2 =*/l  - 2 whence  from  (2.41)  k=X  -X  * . This 

establishes  the  left  decomposition  F = K2Q2  • The  right  decomposition 

T 

F = QjKj  can  be  similarly  established  by  considering  F F in  place  of 
T 

FF1  . 

The  second  part  of  the  result  is  easily  proved  as  follows.  Sup- 
pose now  that  associated  with  the  given  tensor  F there  exists  some 
proper  orthogonal  tensor  Q2  , some  tensor  with  unit  determinant 
and  some  rectangular  cartesian  frame  such  that 


for  some  real  number  k.  Note  that  the  tensors  Q,  , K,  and  the  frame 


X2  are  not  required  to  be  the  particular  ones  used  in  the  preceeding 
proof.  Since  Q2  is  orthogonal,  it  follows  from  (2.50)  that 

EET=!S2kJ  . (2.52) 


T 

whence  in  particular,  the  traces  of  the  two-dimensional  tensors  FF 
T 

and  ^2~2  are  e<lual*  virtue  of  (2.51)  we  now  have  that  necessarily 
1 = FopFap  = 2 + k2  • whence 

k = ±^fT2  . 

The  corresponding  result  for  the  decomposition  F=QjKj  may  be  simi- 
larly established. 

Given  any  plane  volume  preserving  deformation  with  deformation 
gradient  F(x),  we  refer  to  k(x)  defined  by  (2.41)  as  the  associated 
local  amount  of  shear.  Therefore  any  arbitrary  plane  deformation  of 
an  incompressible  material  can  be  viewed  locally  as  a simple  shear  in 
a suitable  direction  with  local  amount  of  shear  k(x)  , followed  or  pre- 
ceded by  a suitable  rotation. 

I s 
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3 • 1 Ellipticity  of  the  Plane  Displacement  Equations  of  Equilibrium 
We  now  introduce  the  relevant  notion  of  ellipticity  without  re- 
stricting ourselves  to  isotropic  materials. 

Consider  a cylindrical  surface  S with  generators  parallel  to 
those  of  the  undeformed  body  and  lying  wholly  within  ft  . Let  C be 
the  curve  along  which  S intersects  II.  Assume  that  C has  a con- 
tinuous curvature,  and  let  ? be  the  arc  length  on  C . Then  C 
may  be  described  by  the  non- singular  parameterization 

C:  x = £ (5)  . 
cr  Of 

If  C is  a coordinate  normal  to  C and  N(5)  is  a unit  vector  normal  to 
C in  the  (Xj,  x.,)- plane,  then  near  a fixed  point  P on  C we  have  the 
orthogonal  curvilinear  coordinate  system  (5,  Ch  permitting  us  to 
write 


xa  = V5,  + CrV5)  (3.1) 

for  any  point  (xj,x2)  in  a two-dimensional  neighborhood  of  P.  The 
mapping  (3.1)  is  locally  one  to  one,  so  that  it  has  an  inverse 

5 ~ ((x j i x2 ) * C = 8(x j»x2)  ( (3.2) 

and  f and  g are  twice  continuously  differentiable  in  a neighborhood 
of  P . Note  that  we  may  take 
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N = 

aw 


Zs_ 

Is*  I 


on  C . 


Now  suppose  that  (u^Xj.x^  , p(xj,x2))  is  a solution  of  the 
plane  displacement  equations  of  equilibrium  (2.3)  and  (2.21)  such  that 
is  once  continuously  differentiable  and  twice  piecewise  continuously 
differentiable  on  II , while  p is  continuous  and  piecewise  continuously 
differentiable  on  fl  . We  set 


V5'C)  = WSJ  + CNjfS)  ' VS,  + CN2(5,)  * 

p(5»  C)  = p(^1(5)  + CNj(5)  , £2(S)  + CN2(5))  , 


and  further  suppose  that,  in  fact,  the  second  order  partial  derivatives 
of  u^  are  all  continuous  across  C except  possibly  for  the  normal 
derivative  9 u^/dC  , and  that  the  first  order  partial  derivative  9p/95  is 
continuous  across  C,  while  the  normal  derivative  9p/9£  may  suffer 
a jump  discontinuity. 

Let 


U = 

a 


P 


(3.3) 


where  [h]  denotes  the  jump  of  a function  h across  C.  Then  one  shows 
easily  that 


[■W'-vyv  • 


(3.4) 


— — — 1/2 

where  N=Vg  = (N*N)  N . We  have  by  the  chain  rule  and  (3.2)  that 


1 


See  Section  1 of 


[1]. 
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P 


a 


which  because  of  the  presumed  smoothness  and  (3.3)  leads  to 

[p#er]=qNa  . (3.5) 

Taking  jumps  in  the  first  two  displacement  equations  of  equilibrium 
(2.21),  and  making  use  of  (3.4),  (3.5)  and  the  assumed  smoothness  we 
get 


C<*Py6UyN6N|3  " F(3aNpq  = 0 on  C • (3.6) 

If  for  all  vectors  N and  nonsingular  tensors  F with  unit  deter- 
minant,  we  define  the  matrix  Qap(N,  F)  through 

Q^l5'9"*l©7l  . (3.7) 

then  is  symmetric  by  virtue  of  (2.  11).  Equation  (3.  6)  can  now  be 

written  in  the  form 


~qFpaNp  on  C * 


(3.8) 


We  also  need  the  "jump  equation"  associated  with  the  remaining 
displacement  equation  of  equilibrium  (2.  3).  We  compute8(detF)/8C  to 
find 


8 1 8 at 


ig-) 

8xpJ 


(3.9) 


where  use  has  been  made  of  (3.2),  the  chain  rule  and  a standard  formula 
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for  the  differentiation  of  a determinant.  Taking  jumps  in  (3. 9)  and  making 
use  of  (2.3),  (3.3)  and  the  presumed  smoothness  leads  to 


But  by  (2.3)  the  jump  in  det  F must  vanish,  whence  (3.  10)  simplifies  to 


Fp‘W°  C . (3.11) 

The  system  of  jump  equations  associated  with  the  displacement  equations 
of  equilibrium  are  (3.8)  and  (3.  11),  and  may  be  regarded  as  three  linear 
homogeneous  algebraic  equations  for  the  jumps  Uq  and  q. 

We  say  that  the  system  of  plane  displacement  equations  of  equili- 
brium is  elliptic  at  the  solution  (u^.p)  and  at  the  point  (Xj,x^)  if  and 
only  if,  for  all  vectors  N/0  , the  system  (3.8),  (3.11)  has  only  the 

trivial  solution  U =0  , q = 0. 

cr 

Consequently  if  the  system  is  elliptic,  the  displacement  field 
uq  will  in  fact  be  twice  continuously  differentiable  at  the  point  under 
consideration  and  the  pressure  p will  be  continuously  differentiable 
there.  If  on  the  other  hand  there  exists  a non-trivial  solution  of  (3.8), 
(3.11)  for  some  vector  N , then  N is  normal  to  a characteristic  curve 

^ /W 

in  the  undeformed  configuration.  These  characteristic  curves  are  the 
only  possible  carriers  of  discontinuities  of  the  kind  admitted  here  in 
u^  and  p , and  ellipticity  precludes  the  existence  of  real  characteristics. 
If  we  set 


m =f:1Na 
a ()<y  p 


(3.12) 


we  can  write  the  system  of  jump  equations  (3.8)  and  (3.  11)  as 
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Q11  Q12  ""M1 

'U1 

Q21  Q22  ”m2 

U2 

■m!  m2  °- 

. q. 

This  system  of  linear  homogeneous  algebraic  equations  for  and  q 
has  only  the  trivial  solution  if  and  only  if 
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the  displacement  equations  of  equilibrium  to  be  elliptic  at  a solution 
(u  p)  and  at  some  point  (x,,x_)  is 

9 1 C 

U,PuFY»F«WQY«<S'S,R«V°  ■ <3-17> 

for  every  vector  N^O  . Finally,  because  of  (3.7)  it  is  clear  that  (3.  17) 
is  equivalent  to 

eo\epuFYXF6uQY6(~'  ~)NaNp^°  for  all  unit  vectors  N . (3.18) 


3 . 2 Specialization  to  Isotropic  Materials 

When  the  material  at  hand  is  isotropic,  we  can  use  (2.4),  (2.  11), 
(2.32)  and  (3.7)  to  simplify  the  necessary  and  sufficient  condition  for 
ellipticity  (3.  18),  which  then  gives 


(‘«‘ftJC»pNlNU>W'«I>  + 2le(,lCc,pNpNx)2'W'(I)/0  . (3.19) 

for  every  unit  vector  N.  Now  let  the  frame  be  principal  for  C,  so  that 


and  evaluate  (3.  19)  in  this  frame.  We  then  find 

(\  \ + X *N*)W'(I)  + 2(X  \ - X *)2N*N*W*(I)^  0 (3.20) 

for  all  unit  vectors  N , as  being  necessary  and  sufficient  for  ellipticity. 
We  will  now  show  that  the  plane  displacement  equations  of 
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equilibrium  are  elliptic  at  a solution  (uQ,p)  and  at  a point  (Xj,X2>  if 
and  only  if 


W'(I)^0  , (I  - 2)  + 1 > 0 , 


(3.21) 


at  the  point  under  consideration;  i.e.  that  (3.21)  is  equivalent  to  (3.20). 


To  show  this,  we  observe  that  since  N is  a unit  vector, 

\ jN2  + \2Ni  = (^  iN2+^2N1)(N1+N2)=X  lN2+X2Nt  + (X  1 +X2)N1N2  ' (3l22) 


so  that  (3. 20)  may  be  written  as 


o 2W'(I)}n|  + [X  j W'(I)  ]N*  + f (X  2 +X  2)W'(I)  + 2(X  2 - X 2 )*W'(I)}n2N2  / 0 (3.23) 


for  all  unit  vectors  N . If  we  set 


E11  = X2W'(I)  ' E22=X1W'(I)  » 


2 ..  2 


E21~E12 


(X  , +x  ?)  _ _ _ _ 

= — L2-£-W/(D  + (Xj  -X‘)V(I),  za=N2  , 


we  can  replace  (3.23)  by 

for  all  z^O,  z a 0 . 
op  o p ~ ~ o 


(3.24) 


(3.25) 


It  has  been  shown  in  Section  2 of  reference  [2]  that  (3.25)  holds  if  and 
only  if 
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E11E22>0 


and 


€ E 


12 


^E11E22 


:>-  1 


(3.26) 


(3.27) 


where 

e = sgn  En=sgnE22  . (3.28) 

Substituting  from  (3.24)  into  (3.26)  we  get  X 2X  2 f W'(I)  }2> 0 
which,  because  X^>0  , is  equivalent  to 

WWO  . (3.29) 

Using  (3.24)  and  (3.28)  in  (3.27)  leads  to 

which  because  of  (2.24)  and  (2.27)  may  in  turn  be  written  as 

2(I"2,lRf)  + 1>0  • (3.30) 

Equations  (3.29)  and  (3.30)  are  what  we  set  out  to  establish. 

A physical  interpretation  of  the  ellipticity  condition  (3.21)  may 

■ ■ 

be  obtained  in  terms  of  the  concept  of  the  local  amount  of  shear  intro- 
duced  in  Section  2.2.  Consider  an  isotropic,  incompressible,  homo- 
geneous, elastic  solid  which  has  a positive  shear  modulus:* 


See  (2.40).  A similar  interpretation  can  clearly  be  given  in  terms  of 
the  local  amount  of  shear  even  in  the  unrealistic  case  when  (2.40)  does 
not  hold. 


(3.31) 


W/(I)>0  for  1*2 


The  first  of  (3.21)  is  now  trivially  satisfied.  If  we  define  the 


function  t by 


T(k)  = 2kW'(2  + k ) , |k|<oo  , 


(3.32) 


then  T(k)  is  easily  shown  to  be  the  shear  stress  corresponding  to  an 
amount  of  shear  k in  a simple  shear  deformation.  The  graph  of  T(k) 
vs.  k described  by  (3.32)  will  be  called  the  response  curve  in  simple 
shear.  Differentiating  (3.32)  with  respect  to  k and  observing  that 
(3.31)  holds  leads  to 


t '(k)  = 2W'(2  + k2)  J 2k2  ■V^2— + 1 
| W'(2  + k ) 


We  therefore  find  that  (3.21)  is  equivalent  to 


'(k)>0  for  k = Jl  - 2 , 


(3.33) 


from  which  we  conclude  that  for  an  isotropic,  incompressible  elastic 
solid  having  a positive  shear  modulus,  the  plane  displacement  equations 
of  equilibrium  are  elliptic  at  a solution  (u^,p)  and  a point  (x^x^)  if 
and  only  if  the  slope  of  the  response  curve  in  simple  shear  at  an  amount 
of  shear  equal  to  the  local  amount  of  shear  is  positive. 

Suppose  for  example  that  the  response  of  a particular  homo- 
geneous, isotropic,  incompressible  elastic  solid  in  simple  shear  is  as 
described  by  Fig.  2.  Then  in  any  plane  deformation  the  displacement 
equations  of  equilibrium  are  elliptic  at  some  point  (Xj,x,)  and  some 
solution  'f  and  only  if  the  local  amount  of  shear  at  that  point  k(Xj,X2>, 
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defined  by  (2.41),  is  such  that  -kQ<  k(Xj,  x2)<  kQ  . 

It  is  apparent  from  the  above  discussion  that  a loss  of  ellipticity 
foi  materials  of  the  type  being  considered  is  dependent  upon  a lo88  of 

inyertibility  of  the  shear  stress  - amount  of  shear  relation  in  simple 
shear. 

Finally,  we  note  from  (3.21)  that  the  undeformed  state  is  ellip- 
tic if  and  only  if  the  infinitesimal  shear  modulus  j3  = 2W'(2)/  0.  This  is 
precisely  the  condition  for  ellipticity  of  the  linearized  displacement 

equations  cf  equilibrium  for  a homogeneous,  isotropic,  incompressible, 
elastic  material. 

3.3  Characteristic  Curves 

If  the  ellipticity  condition  (3.21)  is  violated,  it  follows  that 
there  exists  a unit  vector  N such  that  equality  holds  in  (3.  20).  N will 
then  be  normal  to  a (material)  characteristic,  and  we  now  determine  the 
number  of  possible  characteristics  and  their  inclinations.  To  this  end, 


Nj=  - sin  6 , N?  = cos  6 , 


(3.34) 


so  that  0 is  the  local  inclination  of  the  material  characteristic  to  the 
X j -principal  axis  of  C.  Substituting  this  in  (3.20),  with  equality  hold- 


ing now,  we  find 


(X  jcos20  + X 2sin20)W'(I)  + 2(x 2 - X 2)2sin2ecos2e  W' 


’(I)  = 0 


(3.35) 


We  seek  solutions  0 of  this  equation  in  the  interval  (-f,  f]  . 

Let  us  assume  that  the  infinitesimal  shear  modulus  of  the  ma- 


terial is  positive: 


fl  = 2W'(2)>0  . (3.36) 

We  see  immediately  from  (3.21)  that,  if  the  point  under  consideration 
is  locally  undeformed  (1=2)  in  the  given  deformation,  then  the  dis- 
placement equations  of  equilibrium  are  elliptic  there.  Consequently 
we  need  only  consider  I>2  in  our  search  for  characteristics. 

Suppose  first  that  ellipticity  is  lost  by  virtue  of  the  fact  that  the 
first  of  the  ellipticity  conditions  (3.21)  is  violated.  * Then 

W'(I)  = 0 (3.37) 

at  the  point  (xj,X2>  of  interest  at  the  given  deformation.  We  then  find 
from  (3.35)  that  either  W"(I)  = 0 or  6 = 0,^.  Using  (2.  41 ) and  (3 . 32), 
we  may  state  this  result  as  follows.  Let  k be  the  local  amount  of 
shear.  Then  if  T(k)  = 0,  the  displacement  equations  of  equilibrium  are 
not  elliptic  for  the  given  deformation  at  the  point  under  consideration. 
Furthermore,  we  then  have  two  (material)  characteristics  inclined  at 
angles  0 and  | to  a principal  axis  of  C , except  in  the  particular 
case  when  t '(k)  = 0 as  well,  in  which  case  any  number  of  arbitrarily 
inclined  characteristics  may  exist  locally. 

Now  suppose  that  W'(I)^0  at  the  point  of  interest  and  that  ellip- 
ticity is  lost  by  virtue  of  the  fact  that  the  second  of  (3.21)  has  been 
violated.  Then 

Zvffifj  (I  ~ 2>  + 1 * 0 • (3-38) 

Equation  (3.35)  can  now  be  rearranged  into  the  form  of  a quadratic 


*Note  from  (2.40)  that  this  possibility  does  not  exist  if  the  material  has 
a positive  shear  modulus. 
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equation  for  cos  20  . 


2 2 2 
(a  : -a  ;r 


(X? -x?> 


2 W^I) 


cos  20 


-1J2  J^rx_z\  waal-o 

2 2 Wffi  ~ ' 


(3.39) 


Formally  we  can  write  the  solution  of  this  after  making  use  of  (2.24) 
and  (2. 27)  as 


cos  20  = ■ 


)/w'(i)+i}  {2(i  + 2)w»m/w/m  + 

2(I2-  4)fcw"(I)/W'(I) 


(3.40) 


where  with  no  loss  of  generality  we  have  assumed  that  X j>X^  • 

If  (3.38)  holds  with  equality  (so  that  r(k)^0  , T#(k)  = 0 at  the 
local  amount  of  shear  k)  we  find  two  values  of  0 in  the  interval 

from  (3.40),  whence  two  characteristics  exist.  Equation  (3 . 40) 


now  simplifies  to 


“•  26  * - ,/m  • 


(3.41) 


which  because  of  (2.24)  and  (2.27)  (and  since  Xj>X2)  leads  to 


cos  20 


(3.42) 


whence 


tan  0 = ±X 


(3.43) 


Suppose  the  spatial  characteristic  corresponding  to  this  material  char- 
acteristic is  inclined  at  an  angle  a to  the  X j -principal  axis  of  G . 

It  can  be  shown  that 


^ 2 

tana  = — tan  9 , (3.44) 

A 1 

so  that  (2.24),  (3.43)  and  (3.44)  give 

tano  = ±p  . (3.45) 

A 1 

Because  of  (2.33),  a is  also  the  inclination  to  the  corresponding  prin- 
cipal axis  of  the  Cauchy  stress  tensor. 

If  however,  strict  inequality  holds  in  (3.38)  (so  that  T(k)j*0  , 
kT(k)r/(k)<0  at  the  local  amount  of  shear  k)  (3.40)  gives  us  four  values 
of  8 which  in  turn  implies  the  existence  of  two  pairs  of  characteristics. 
Clearly,  each  pair  is  positioned  symmetrically  with  respect  to  the  prin- 
cipal axes  of  C.  In  what  follows  we  will  have  need  for  the  inclinations 
a of  the  corresponding  spatial  characteristics  to  the  X j -principal  axis 
of  G (Xj>X2).  From  (2.  24),  (2.27),  (3. 40)  and  (3. 44)  we  have 

COS  2a  = - 2(1  - 2)W»(I)/ W'(I)  +1 } f 2(1  4 2)W«(I)/ W'(I)  4 U )*  _ 

2 (I2-  4)fc\V"(I)/W'(I) 


(3.46) 


4.  1 Weak  Formulation  of  Problem 


In  the  derivation  of  the  classical  field  equations  of  elasticity  the 
displacement  field  ju  and  stress  field  £ are  assumed  to  satisfy  certain 
smoothness  requirements.  There  are,  however,  some  physical  prob- 
lems in  which  these  conditions  are  not  met,  so  that  in  order  to  study 
them  one  would  be  forced  to  relax  the  smoothness  demanded  of  the  field 
quantities.  It  may,  for  example,  be  necessary  to  require  only  that  the 
displacement  field  £(x)  be  continuous  and  piecewise  continuously  dif- 
ferentiable on  ft  , while  the  nominal  stress  field  £(x)  and  the  pressure 
field  p(x)  are  to  be  piecewise  continuous1  on  ft  . Clearly,  the  global 
balance  laws  continue  to  be  meaningful  even  under  these  smoothness 
conditions,  but  one  must  re-examine  the  validity  of  the  local  field 
equations. 

Of  particular  physical  interest  is  the  case  wherein  the  field 
quantities  possess  the  classical  degree  of  smoothness*”  everywhere 
except  on  one  or  more  regular  surfaces  within  the  body.  This  would, 
for  example,  describe  an  idealized  model  for  shear  bands.  To  formu- 
late this  problem,  we  suppose  that  there  is  a surface  S in  ft  such 
that  £ , F and  p are  continuously  differentiable  everywhere  in  ft 
except  on  S , and  such  that  £ , F and  p suffer  finite  jump  disconti- 
nuities across  it.  The  displacement  u(x)  is  presumed  to  be  continuous 

* We  return  momentarily  to  the  three-dimensional  case  in  this  section. 

2 

See  Section  2.1. 
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everywhere  in  Si  . The  possibility  of  the  breakdown  of  ellipticity  of  the 
governing  equations  suggests  that  solutions  of  this  type  to  the  equations 
of  finite  elastostatics  may  emerge  in  some  circumstances. 

Going  through  the  usual  arguments,  * one  finds  from  the  global 
equilibrium  of  forces  that 

divo  = 0 on  ft  - S (4.1) 

<■>*/ 

and 

[o]  + N = 0 on  S , (4.2) 

while  from  the  global  equilibrium  of  moments  we  have 

0FT=  F0T  on  ft  - S (4.3) 

and 

y(x)x[o]  + N = 0 on  S . (4.4) 

Equation  (4.2)  says  that  the  nominal  tractions  are  continuous  across 
r ,+  + + 

S . Here  [aj  = 0-0  where  a and  a are  the  limiting  values  of  a 
(presumed  to  exist)  as  a pointon  S is  approached  from  each  side,  and  N 
is  is  a unit  normal  to  S . Equations  (4.  2)  and  (4.  4)  are  referred  to  as 
jump  conditions.  Note  that  (4.4)  is  trivially  satisfied  once  (4.2)  is. 
Incompressibility  likewise  leads  to 

detF  =1  on  ft  - S . (4.5) 


See  Chadwick 


[7], 


page  114. 


„<■■■.  i. 
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Such  a surface  S carrying  jump  discontinuities  in  F,  O and 
p which  conform  with  the  jump  condition  (4.2),  while  maintaining  con- 
tinuous displacements  across  it  is  called  an  "equilibrium  shock",  or 
an  "elastostatic  shock"  in  the  particular  case  when  the  body  is  com- 
posed of  an  elastic  material. 

4. 2 Piecewise  Homogeneous  Elastostatic  Shocks 

To  investigate  many  of  the  local  issues  related  to  elastostatic 
shocks,  it  is  sufficient  to  consider  the  case  in  which  S is  a plane  and 
the  deformation  gradient  F is  constant  on  either  side  of  S . From 
here  on  we  shall  be  concerned  with  such  a situation  within  the  context 
of  plane  deformations*  of  an  incompressible  elastic  solid,  so  that  we 
may  take  S to  be  a plane  parallel  to  the  generators  of  the  body. 

The  corresponding  problem  for  a compressible  elastic  solid 
was  investigated  by  Knowles  and  Sternberg  [3],  In  this  section,  we 
formulate  the  problem  governing  the  existence  of  an  elastostatic  shock 
in  the  incompressible  case  in  a manner  entirely  analogous  to  [ 3 J . 

Suppose  that  the  middle  cross-section  of  the  body  we  are  deal- 
ing with  occupies  the  entire  (x x^J-plane  II  in  its  undeformed  configu- 
ration. Let  X be  a fixed  rectangular  cartesian  coordinate  frame  and 
let  £ be  the  straight  line  through  the  origin  of  X with  unit  direction 

vector  L . Thus 
*** 

£ i kq=  Lq5  , -co<5<oo  . (4.6) 


We  leave  the  three-dimensional  introduction  to  elastostatic  shocks  of 
the  last  section  and  return  to  plane  deformations  from  here  on. 


Let  N be  the  unit  vector  normal  to  J.  obtained  by  a counterclockwise 

+ 

rotation  of  L . Let  FI  and  II  be  the  two  open  half  planes  into  which 

+ 

^ divides  II  , with  II  being  the  one  into  which  N points.  (See  Fig.  1.  ) 
Now  consider  the  piecewise  homogeneous  plane  deformation 


where  F and  F are  constant  tensors  such  that 


The  nominal  stresses  associated  with  the  deformation  (4.7)  are 


aw(F)  + + 

ap-^-pF, 


8W(F) 


Clearly,  the  equilibrium  equations  (4.  1)  are  satisfied  if  and  only  if  p 
and  p are  constants. 

If  we  are  to  view  the  line  £ as  the  intersection  of  an  equilibrium 
shock  S with  the  cross-section  II  , then  according  to  Section  4.  1 we 
need  to  impose  displacement  and  traction  continuity  requirements  across 
. Because  of  (4.7)  the  requirement  of  a continuous  displacement  field 
is  equivalent  to 
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which  in  view  ol  (4.6)  reduces  to 


(.  La  = F A La  . 

aP  P orP  P 


By  (4.9),  we  have  traction  continuity  (4.2)  if  and  only  if 


8W(F)  + + _! 

- pF 


8F 


ap 


pa 


> v< 


9W(F) 


> N 


P ‘ 


(4. 11) 


(4.  12) 


If  the  deformation  field  (4.7),  subject  to  (4.8),  together  with 

+ - + - 
real  constants  p and  p conform  with  (4.  1 1)  and  (4.  12),  and  if  F^  F , 

then  we  refer  to  the  corresponding  elastostatic  field  as  a piecewise 
homogeneous  elastostatic  shock.^  The  line  iC  will  be  referred  to  as 
the  material  shock-line.  Figure  1(b)  displays  the  images  of  the  three 
rectangles  shown  in  Fig.  1(a)  under  a typical  mapping  (4.7)  in  the  pres- 
ence of  such  a shock. 

In  order  to  examine  questions  related  to  the  existence  of  piece- 
wise  homogeneous  elastostatic  shocks  we  pose  the  following  problem. 

+ + X 

Given  a constant  tensor  F with  det  F = 1 and  a real  constant  p , de- 

termine  a constant  tensor  F with  det  F = 1 (F/ F)  and  a real  constant 

p such  that  (4.  11)  and  (4.  12)  hold. 

Equation  (4.  11)  may  be  solved  as  follows.  Let  , which  we 

shall  refer  to  as  the  spatial  shock-line,  be  the  image  of  / under  the 

mapping  (4.7).  Let  and  be  the  two  half  planes  into  which  II 

and  IT  map  by  virtue  of  (4.7).  Suppose  is  the  unit  direction  vector 

of  such  that  the  unit  normal  n to  obtained  by  rotating  / 

+ 

counterclockwise  points  into  W*  . (See  Fig.  1.)  Without  any  loss  of 

*Note  from  (4.  12)  that  if  p^p  then  necessarily  F^  F . 


J 
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generality  the  inclinations  * and  0 of  the  shock-lines  and  J. 
relative  to  the  x^-axis  may  be  confined  to  the  interval  f - ^ ^ j . 

One  can  show  readily  that,  given  a constant  tensor  F with 

+ 

det  F = 1 , (4.  11)  will  hold  for  a tensor  F with  unit  determinant  if  and 

<v  fV 

only  if 


+ 

F ft=  (6  n )F 

op  cry  a y yP 


(4.  13) 


for  some  real  number  k . We  omit  the  derivation  of  this  result  as  it 
parallels  exactly  the  corresponding  derivation  in  the  compressible  case 
contained  in  [3],  Let  X'  be  the  rectangular  cartesian  frame  obtained 
by  rotating  the  frame  X counterclockwise  through  an  angle  0 . The 
base  vectors  associated  with  X'  are  then  £ and  n . Expressing  (4.13) 
in  the  frame  X'  we  have 


;X' 

12 


FX' 
* 22 


-1 

r i 

r 

i 

H 

0 

1 

J 

L J 

L 

+X' 

11 

+x' 

21 


+ 

F 


X ' 
12 


(4.  14) 


Accordingly,  the  deformation  on  II  may  be  viewed  as  being  equivalent 

+ 

to  the  deformation  on  II  followed  by  a simple  shear  parallel  to 
with  an  amount  of  shear  k . 

We  may  now  pose  the  following  problem  which  is  equivalent  to 

+ 

the  one  posed  earlier.  Given  a constant  tensor  F with  unit  determi- 
nant and  a real  constant  p , determine  real  numbers  p , k(^0)  and 
0 such  that  (4.  12)  holds  with  F defined  by  (4.  13).  Here 

we  have  omitted  k = 0 since,  by  (4.  13)  we  see  that  this  corresponds  to 
the  shockless  state  F = ^ . 


Finally,  we  note  that  since  the  traction  continuity  condition  (4.  12) 
imposes  only  two  scalar  restrictions  on  the  three  parameters  0 , k 
and  p , one  would  anticipate  that  if  there  exists  an  elastostatic  shock 
corresponding  to  a given  F and  p , then  in  fact  there  exists  a one  - 
parameter  family  of  shocks. 
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5.  1 Weak  Piecewise  Homogeneous  Elaatostatic  Shocks 

We  now  specialize  the  problem  posed  in  the  general  setting  of 

Section  4.2  to  the  first  of  two  simpler  cases.  Here  we  confine  attention 

to  elastostatic  shocks  that  are  weak  in  the  sense  that  the  departure  of 
- + 

F from  F is  small.  Motivated  by  the  remarks  at  the  end  of  the  pre- 
vious section,  we  assume  here  that  there  exists  a one-parameter  family 
of  shocks,  corresponding  to  the  given  F and  p , depending  on  the 
parameter  h and  sufficiently  smooth  near  k = 0 . Specifically,  we 
suppose  that  there  are  functions  0(x)  , p(x)  both  twice  continuously 
differentiable  in  a neighborhood  of  k = 0 , such  that  F defined  by  (4.  13) 
together  with  p(x)  conforms  with  the  traction  continuity  requirement 
(4.  12),  Since  from  (4.  13)  we  have  that  F = F when  h = 0 , we  may  use 
h as  a measure  of  the  departure  of  F from  F . Accordingly  k will 
be  referred  to  as  the  shock-strength  parameter. 

We  first  record  the  following  kinematic  results  which  are  estab- 
lished in  [3  ] . Let 


c=  FL  = FL  . 


Then 


*ot=  cFa(3Lp=  cFa|3Lp  ' 


no=cFpuNp=c^paNp  * 


(5.  1) 

(5.2) 

(5.3) 
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Vp  +nanp  = 6ap  ' La"eapNp  * . (5-4) 

If  0(h)  and  p(h)  exist  as  described  above,  it  follows  that  L , N . f , 

• ♦ 
n , F and  c are  all  h dependent  whence  we  write  L(k)  , N(h)  , f(n)  , 

n(n)  , F(h)  and  c(k). 

Because  of  the  presumed  smoothness  of  0(h)  we  have  the  fol- 
lowing Taylor  expansions,  where  a prime  denotes  differentiation  with 
respect  to  h , 


i(H)  = i(0)  + i'(0)H+o(H)  , £(*}  = n(0)+n'(0)H+o(K)  , (5.  5) 1 

N(h)=N(0)  + N'(0)h+o(h)  . 

Equation  (4.  13)  now  gives 


wk)=wkvo)\(o)^y0+o(h)  • 

(5.6) 

Fjp1' <*>  = -■‘\(0>np(0)F<t-y1  + o(k)  , 

(5.7) 

where  we  have  also  used  (5.4).  This  enables  us  to  write  the  following 
Taylor  expansion 

8W(F(x|)  8W(F)  82W(F)  + 

The  Taylor  expansion  of  p(n)  leads  to 
^Whenever  we  write  o(h)  , we  mean  o(h)  as  h-*0  . 
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p(k)=P(0)  + kp/(0)  + o(k) 


(5.9) 


Using  (5.6),  (5.7)  and  (5.9)  and  evaluating  the  traction  continuity  con- 
dition (4.  12)  to  leading  order,  gives 


p = p(0)  . 


(5.  10) 


Consequently  we  may  write  (5.9)  as 


p(x)  = p + np  '(0)  + o(x) 


(5. 11) 


We  now  return  to  the  traction  continuity  condition  (4.  12)  and  re-evaluate 
it  to  leading  order  using  (2.11),  (5.5),  (5.7),  (5. 8)  ana  (5.  11).  This 


leads  to 


4y°)no<°>*pY‘K<0> ■ ° - (5-l2) 

which  are  two  scalar  equations  for  p'(0)  and  0(0)  . 

5.2  A Necessary  Condition  for  the  Existence  of  a Weak  Shock 

We  now  derive  a necessary  condition  for  equation  (5.  12)  to  have 
a solution  p'(0)  , 0(0)  . We  have  from  a Taylor  expansion  of  (5.3)  that 


(5.  13) 


Equations  (5.  12)  and  (5.  13)  lead  to 


(F)lv(°)nw(0)nv(0>^6Fvp”  P '(0)no(0)  + J I Y(0)nY(0)no(0)  = 0 . (5.14) 


"0fPy6~  Y 


But  since  i is  perpendicular  to  n we  have  fo(0)na(0)  = 0 , whence 


(5.  14)  simplifies  to 
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c^yt(E»,v<0»"„<0>''v«»F„6Fv(3.p'(0)„a(0)  = 0 . (5. 

Multiplying  (5.  15)  by  na(0)  and  maxing  use  of  the  fact  that  n is  a unit 
vector  leads  to 


f'(0,*C«M(£>\<0K.(0V0>V0,^,A|l  . (5. 

Alternatively,  multiplying  (5.  15)  by  f^(0)  gives 

%Pv6^(0V0>V0V0>^vp  = ° , ,5. 

by  virtue  of  the  fact  that  f- n = 0 . Using  (5.3)  and  (5.4)  in  (5.  17)  leads 
to 


*Trve\MFy7rFoACQ-pv6(E)Np(0)N6(0)Nv(0)Nu(0)  = 0 , (5.18) 

which  because  of  (3.7)  can  be  equivalently  written  as 

*o^c^Fv\F6(aQy6(^(0)'£)No(0)Np(0>  = 0 • (5.19) 

Equation  (5.  19)  must  necessarily  hold  if  a one  parameter  family  of 
elastostatic  shocks  of  the  type  being  considered  is  to  exist.  On  com- 
paring with  (3.  18),  we  see  that  (5.  19)  implies  a loss  of  ellipticity  of  the 
displacement  equations  of  equilibrium  on  n at  the  given  F and  d 
We  therefore  have  the  following  results 

Theorem  1.  A necessary  condition  for  the  exist- 
ence of  a one-parameter  family  of  elastostatic 
shocks,  of  the  kind  under  consideration,  is  that 
the  displacement  equations  of  equilibrium  suffer 
a loss  of  ellipticity  at  the  given  deformation  and 
pressure  on  n . Furthermore,  in  the  weak  shock 
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limit  (h-»0)  the  material  shock-line  and  the 
spatial  shock-line  tend  respectively  to  a ma- 
terial and  spatial  characteristic  associated 
with  n . 

The  corresponding  result  was  obtained  by  Knowles  and 

Sternberg  [3]  in  the  case  of  compressible  elastic  materials. 

In  the  event  that,  corresponding  to  a given  F and  p a one- 

parameter  family  of  shocks  of  the  type  being  considered  exists,  the 

jumps  of  various  physical  quantities  across  the  shock  can  be  easily 

+ + 

determined  to  leading  order  in  terms  of  the  given  F , p and  the  pre- 
sumably determinable  (from  (5.12))  0(0)  , p'(0)  . We  now  determine 
some  of  these  jumps. 

(i)  The  jump  in  energy  density  [w]  + 

The  Taylor  expansion  of  W(F(k))  about  h = 0 , together  with 
(2.9)  and  (5.6)  leads  to 

W(F(K))  = W(F)+K{;J:apnp(0)fcr(0)  + pia(0)na(0)}  + o(K)  . (5.20) 

Since  is  perpendicular  to  n we  can  drop  the  last  term  in  (5.20)  to 
get 

[W];  = Hta(0)fa(°)  + o(K)  , (5.21) 

where  we  have  set 

t(n)  = Tn(x)  . (5.22) 

As  a consequence  of  (2.6),  (4.2),  (5.3)  and  displacement  continuity,  we 
see  immediately  that  t(H)  = Tn(x)  = Tn(x)  which  implies  the  continuity  of 
the  Cauchy  traction  vector  across  . 


(ii)  The  stress  jumps  t TCrp 3 + 


From  (2.9)  we  have  that 


Tcrp'aFap  Fpv‘p6c*P  * 


(5.23) 


+ 9W(F)  + 


ToP=  9Fop  FPv'PW  » 


(5.24) 


which  together  with  (2.  11),  (5.6),  (5.8),  (5.  1 1 ) and  (5.  22)  leads  to 


[T«PI;-{‘;mJ,6'El^6^(3u'v<0)n.(0)-P'(0)6!,|3+Pna(0)'p<0» 


+ ta<O)f0(O)}  + o(K)  . 


(5.25) 


(iii)  The  jump  in  the  normal  stress  acting  on  a plane  perpendi- 


cular  to  ^ 


Consider  the  plane  perpendicular  to  the  spatial  shock  so  that  the 


normal  to  this  plane  is  £ . The  jump  in  the  normal  stress  acting  on 


this  plane  across  the  shock-line,  [tjjJ  * is 


r X'i+  + . . - , . 

tTlll-  = T«pVp"  T«pV< 


op  Or  p crP  or  P ' 


(5.26) 


which  because  of  (5.25)  and  the  perpendicularity  of  the  vectors  t , and 


n can  be  written  to  leading  order  as 


fjV  = -«{ to(°)'a(0)  - p'(°) 


tcauv6(E)^5^Pu'v(O)V0)'=.(0)nT.<O,)  + '’(K)  ' <5'27) 


In  view  of  (5.4)  and  (5.  1 , this  leads  to 


-45- 


6.  1 Finite  Elastostatic  Shocks  in  Isotropic  Incompressible  Materials 
We  now  return  to  shocks  of  finite  strength,  but  assume  the  ma- 
terial at  hand  to  be  isotropic.  Substituting  (2.31)  in  (4.  12)  and  making 
use  of  (2.4)  and  (5.3)  leads  to 

= {2W'(f)Ga(J-  . (6.1) 

where 


1 = FcrpFcrp  = GQfCr  ' * = FapFop  = Gocr  • 


(6.2) 


Clearly,  (6. 1)  is  simply  a statement  of  the  fact  that  the  Cauchy  traction 
is  continuous  across  the  spatial  shock.  The  original  problem  concern- 
ing the  existence  of  elastostatic  shocks  can  now  be  posed  as  follows: 

given  a constant  tensor  F with  unit  determinant  and  a real  constant 
+ - / 

p , determine  real  numbers  p , k(7=0)  and  0 such  that  (6.  1)  holds 
with  G given  by  (4.13),  (6.2). 

If  we  express  (6.1)  in  terms  of  its  components  in  the  frame  X', 
we  have 

2G^2V(I)  = 2G^2V(i)  , (6.3) 

2G22  W'(I)  - p = 2G*  V(i)  - p . (6.4) 

As  observed  earlier,  (6.3)  and  (6.4)  impose  only  two  scalar  restrictions 
on  the  three  quantities  0 , k and  p . Furthermore  since  p enters 
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only  in  (6.4),  and  there  too  only  linearly,  we  may  consider  (6.3)  and 
(6.4)  separately  i.e.  if  there  are  numbers  h and  0 such  that  (6.3) 
holds,  then  there  certainly  is  a third  number  p such  that  (6.4)  holds 
as  well.  The  existence  of  an  elastostatic  shock  therefore  depends  on 
whether  there  are  numbers  k and  0 such  that  (6.3)  holds. 

To  pursue  this  question  further,  we  need  the  components  of 

+ 

G and  G in  the  frame  X*  . With  no  loss  of  generality  let  us  take  X 

+ 

to  be  a principal  frame  for  G . Then 

0 


, X 2 - 1 


(6.5) 


By  the  change  of  frame  formula  for  tensors  we  have 


<*XU3X-T 


- R , 


(6.6) 


R = 


cos  0 sin0 


-sin0  co 8 0 


(6.7) 


Equations  (6.5)  - (6.7)  lead  to 


+ v' 
GX  = 


2 2 2.2 
X i +x2  x i -x  2 
= + 5 cos  20 


2 2 
*1**2 


sin  20 


sin  20 


x2+x2 
K 1 *2 

— T~ 


x2-x2 

^1  K2 


cos  20 


(6.8) 


2 
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If  we  now  set 


2 2 
X1'X2 
P=4 


XI+A2 


(6.9) 


we  can  write  (6.  8)  because  of  the  second  of  (6.  5)  as 


+ v ' 
G X = 


J i-p2 


1 + P cos  20  - (3  sin  20 


-|3sin20  1 - p cos  20 


(6.  10) 


It  is  clear  from  (6.  5)  and  (6.9),  that  the  value  of  p alone  suf- 


+ X 


fices  to  determine  G completely,  and  in  this  sense  p is  a measure 


of  the  deformation  on  IT  . Note  that  because  X^>0  , (6.9)  implies  that 


1>P>-1  . 


(6.  11) 


Furthermore,  we  have  P = 0 if  and  only  if  the  part  of  the  body  occupy - 
+ 


ing  II  in  its  undeformed  configuration  remains  undeformed  under  the 
mapping  (4.7). 

We  now  find  from  (4.  14),  (6.  2)  and  (6.  10)  that 


X'  1 


J 1 - P2 


1 + pcos  20  - 2k  Psin  20 


+ h (1  - Pcos  20) 


- Psin  20 
+ k(1  - pcos  20) 


- Psin  20  + k(1  - Pcos  20)  1 - Pcos  20 


(6.  12) 


is  then  a proper  orthogonal  tensor. 


■ 


i 
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and  from  (6.  2),  (6.  10)  and  (6.  12)  that 


(6. 13) 


(6.  14) 


Returning  to  the  traction  continuity  requirement  (6.3)  with  (6.  10), 
(6.  12)  - (6.  14)  we  find 


We  may  now  pose  the  problem  as  follows:  given  a number  p in  (-1,1), 
find  numbers  k^O  and  0 in  f - ^ ^ J such  that  (6.  15)  holds. 

If,  for  the  given  p in  (-1,1)  and  any  0 in  [ " ^ ^ j there 

does  not  exist  a root  h/0  to  (6.15),  the  material  is  incapable  of  sus- 
taining an  elastostatic  shock  corresponding  to  the  given  deformation 

+ 

associated  with  p on  II  . On  the  other  hand  if,  for  the  given  p in 
(-  1,1)  and  for  some  0 in  [*-^,  there  is  a root  k/0  to  equation 

(6.  15),  then  there  exists  a corresponding  elastostatic  shock.  There- 
fore, we  now  investigate  the  possibility  that  (6.  15)  has  a root  h^O 
for  all  values  of  0 and  P such  that  ‘|S0S|  , |p|<l  . 

Finally  we  observe  from  (6.  15)  that  if  for  some  pair  (0  , p)  , 
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I P I < 1 and  , there  exists  a root  k to  equation  (6.  15),  then 

(i)  -x.  is  a root  of  (6.  15)  for  the  values  (-0  , (3)  , and 

(ii)  - k is  a root  of  (6.  15)  for  the  values  - 0 , - (3)  . 

It  therefore  follows  that,  as  far  as  the  issue  of  existence  is  concerned, 
we  may  in  fact  restrict  0 to  f~0,  and  (3  to  [0,  1).  If  we  define  the 
set  G by 

G ={(0,(3)  1 0*0*!  , 0*|3<1}  , (6.16) 

we  need  to  look  at  the  question  of  the  existence  of  a root  h^O  to 
equation  (6.  15)  for  every  (0,  (3)  in  G . 

6. 2 Some  General  Results 

We  now  establish  some  general  results  concerning  the  existence 
of  elastostatic  shocks,  valid  for  a homogeneous,  isotropic,  incompres- 
sible elastic  solid  which  has  a positive  shear  modulus. 

We  first  make  the  following  preliminary  observation.  If  (3  = 0 
or  0 = 0 or  0=^  , the  only  root  of  (6.  15)  is  k =0  . This  follow? 
directly  from  (6.  15)  because  of  (2.40).  Consequently,  for  a material 

of  the  type  we  are  considering,  no  elastostatic  shock  is  possible  if  the 

+ 

part  of  the  body  occupying  II * is  undeformed,  nor  can  any  spatial 
shock-line  be  inclined  at  0 or  ^ to  the  principal  axes  of  G . We 

O 

may  now  restrict  attention  to  the  interior  G of  the  set  G: 

G ={(0.P)|o<0<!  , 0<  P<  1}  . (6.17) 

If  we  set 

*We  assume  from  here  on  that  (2.40)  holds. 


on  G , (6.19) 

we  can  write  (6.  15)  using  (6.  13),  (6.  18)  and  (6.  19)  as 

bW/(I)  = (b  + ck  )W/(l'+  2 bn  + ck2)  . (6.20) 

Clearly 


b<0 

on 

O 

G , 

(6.21) 

c>0 

on 

0 

G . 

(6.22) 

A 

© 

Choose  and  fix  a point  (0 

. P) 

in  G . 

At  this  fixed  value  of 

0 and  p we  define  the  function  h by 

h(x)  = (b  + cx)W,(I  + 2bx  + cx2)-bW/(I)  for  |x|<oo  , (6.23)1 

where  1 , b and  c are  given  by  (6.  13),  (6.  18)  and  (6.  19)  respectively 
evaluated  at  (0  , p)  . If  the  plane  strain  elastic  potential  W(I)  is  twice 
continuously  differentiable  on  Is  2 , as  we  have  implicit. y assumed,  it 
follows  that  h(x)  is  continuously  differentiable  on  (-  oo  , oo).  If  there 
exists  an  equilibrium  shock  corresponding  to  the  homogeneous  defor- 

a + a 

mation  associated  with  p on  II  and  inclined  at  an  angle  0 to  the 
y^-axis,  it  is  necessary  and  sufficient  that  h(x)  have  a zero  at  some 

*From  (6.  13),  (6.  18),  (6. 19)  and  (6. 22)  we  have  that 

1 + 2bx  + cx2  = 2 + H?-t.5r£)...  + ^-c~— 1—  i 2 for  all  |x  | < co  . 

C C ’ 


c = c(0,P)  =• 


- Pcos  2j 
(1-p2)* 


0 . The  zero  of  h(x)  gives  the  shock  strength  k . 


Because  of  (2.40),  (6.21)  and  (6.22)  we  see  from  (6.23)  that 

h(0)  = 0 , (6.24) 

h(-£)=-bW'(I)>0  , -£>0  . (6.25) 

V.  C 

It  now  follows  from  (6.24),  (6.25)  and  the  smoothness  of  h(x)  that: 

(i)  if  h'(0)<0  , then  there  exists  a zero  of  h(x)  other  than 
x = 0 , 

(ii)  if  there  exists  a zero  other  than  x = 0 of  h(x)  , then  there 

exists  a number  such  that  h,(K:f.)  = 0.  Furthermore, 

■f  O 

since  h'(-  b/c)  = cW'(I-b'/c)>0  we  have 


• (6.26) 

Because  of  the  remarks  made  before  (6.24),  we  may  interpret 
(i)  and  (ii)  as  follows: 

(a)  h/(0)<0  is  sufficient  to  ensure  the  existence  of  an  elasto- 
static  shock  corresponding  to  the  deformation  associated 

A + A 

with  P on  IT  with  spatial  shock  inclination  0 . 

(b)  If  an  elastostatic  shock  as  just  described  is  to  exist,  then  it 
is  necessary  that  h'fa^sO  for  some  number  . 

This  leads  to  the  main  results  of  this  section  which  we  now  es- 
tablish. We  first  introduce  the  following  terminology.  Recall  from 
(2.40)  and  (3.21),  that  a loss  of  ellipticity  of  the  displacement  equations 
of  equilibrium  can  occur  at  some  deformation  and  at  some  point  if  and 
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(I-Z)  + UO  . (6.27) 

If  ellipticity  is  lost  because  in  fact  the  strict  inequality  holds  in  (6.27), 

we  say  that  a strict  failure  of  ellipticity  has  occurred.  Since  four 

characteristic  curves  exist  in  this  case  (see  Section  3.3)  one  may  say 

that  the  displacement  equations  of  equilibrium  are  hyperbolic  at  such 
+ 

a deformation  on  II  . 

Theorem  2.  A strict  failure  of  ellipticity  of  the  dis- 
placement equations  of  equilibrium,  at  the  given  homo- 
geneous deformation  and  pressure  on  II  , is  sufficient 
to  ensure  the  existence  of  a corresponding  elastostatic 
shock  in  a homogeneous,  isotropic,  incompressible, 
elastic  solid  with  a positive  shear  modulus.  * 

Proof: 

By  hypothesis,  the  given  deformation  gradient  F is  such  that 

A 

the  associated  value  of  (3  , say  p , given  by  (2.4),  (6.5)  and  (6.9)  con- 
forms with  the  inequality 

^U-2m<0  . (6.28) 

W'(I> 


where  by  (6.13) 


+ 

I = 


(6.29) 


Note  from  (2.40)  and  (6.28)  that  necessarily  1^2  , whence  P^O  . We 

A 

now  choose  the  value  of  0 as 


This  and  the  following  result  can  be  readily  modified  for  materials  whose 
shear  modulus  is  not  always  positive. 


We  will  show  that  corresponding  to  the  given  homogeneous  deformation 
on  n , such  that  the  associated  value  of  p conforms  with  (6.28)  and 

A 

(6.29) ,  there  exists  an  elastostatic  shock  at  the  inclination  0 giver,  by 

(6.30) . 

According  to  (6.23)  and  statement  (a)  following  (6.26),  we  need 
only  show  that 

2b2W"(I)  + cW,(I)<0  at  ($,(3)  , (6.31) 

in  order  to  establish  this.  Using  (2.40)  we  may  write  (6.28)  alternate- 
ly as 


2b2W*(I)  + cW'(f)  + 


(6.32) 


A A 

where  b and  c are  defined  by  (6.  18),  (6.  19)  and  evaluated  at  (0,  p)  . 
Using  (6.  18),  (6.  19)  and  (6.  29)  in  (6. 32)  we  find 


2b2W*(I)  + cW'(I)  - 


A 2 

-&T 


(l-p  )(I  -2) 


- 1 -J  l-p* 

cos  20 — 


? W'(I)<  0 , (6.33) 


which  because  of  (6.30)  reduces  to  (6.31),  which  in  turn  establishes  our 
result. 

Theorem  3.  A necessary  condition  for  the  existence 
of  a piecewise  homogeneous  elastostatic  shock  in  a 
homogeneous,  isotropic,  incompressible,  elastic 


^Since  P^O  , |p|<l  this  defines  a real  angle  0 in  (0,-^)  . 


W 
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solid  with  positive  shear  modulus,  is  that  the  dis- 
placement equations  of  equilibrium  suffer  a loss 
of  ellipticity  at  some  homogeneous  deformation. 


Proof: 


By  hypothesis  there  is  a point  (0  , (3)  in  G such  that  there 
exists  an  associated  elastostatic  shock.  By  statement  (b)  following 
(6.26)  then,  there  is  a real  number  k*  such  that 


h'(H*>  = 0 


(6.34) 


where  h(x)  is  given  by  (6.23)  with  b , c and  I evaluated  at  (0  , (3) 
Equations  (6.23)  and  (6.34)  give  that 


Let 


cW'(I  + 2 bn  *+  cn  2 ) + 2(b  + ch  *)2  W*U  + 2bn  *+  ck  2)  = 0 . 


I*  = I +2bn*+cn2  , 


(6.35) 


(6.36) 


so  that  we  have  from  (6.3  5)  that 


(I*-  2){cW'(I*)  + 2(b  + cn  *)2W#(I#)}  = 0 . 


It  follows  from  this  that 


(I*-2){cW'(I*)  + 2(b  + CH*)2W^I*)} 


(i-P  ) 


cos  20  - 


1 -S\  - P2 


P 


' 2 

y w'dj  , 


(6.37) 


since  by  virtue  of  (2.40)  and  (6.  11)  the  right  hand  side  of  (6.37)  it 


• ‘ v > !te  sfK 
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non-negative.  Using  (6.13),  (6.18),  (6.19)  and  (6.36)  in  (6.37)  leads  to 

(b  + CK*)2{2W%)(V2)  + W%)}sO  , (6.38) 

which  because  of  (2.40)  and  (6.26)  gives 

frffc  -j  (V  2) +1*0  . (6.39) 

This  implies  a loss  of  ellipticity  of  the  displacement  equations  of  equili- 
brium at  a homogeneous  deformation  in  which  the  deformation  gradient 
F is  such  that  = 1^  . 

To  summarize,  we  have  shown  that  for  the  type  of  material  at 
hand,  a strict  loss  of  ellipticity  at  the  given  deformation  is  sufficient 
to  ensure  the  existence  of  a corresponding  elastostatic  shock.  On  the 
other  hand,  a loss  of  ellipticity  at  some  homogeneous  deformation  is 
necessary,  if  an  elastostatic  shock  is  to  exist. 

We  draw  attention  to  the  fact  that  Theorem  2 does  not  imply 
that  if  ellipticity  is  strictly  lost  at  the  given  deformation  then  the  cor- 
responding configuration  of  the  body  must  involve  a shock.  Rather,  it 
claims  that  such  a configuration  is  available.  There  is  also  a shock- 
less configuration  available  corresponding  to  the  root  h = 0 of  (6.20). 
Likewise,  a loss  of  ellipticity  at  the  given  deformation  is  not  necessary 
for  a corresponding  elastostatic  shock  to  exist.  In  a boundary-value 
problem  that  we  have  studied,  the  results  of  which  will  be  reported  in  a 
separate  paper,  we  encountered  configurations  of  a body  involving  elasto- 
static shocks  such  that  the  displacement  equations  of  equilibrium  were 
elliptic  on  both  sides  of  the  shock-line. 
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i ' 

7.  1 Dissipativity  Inequality 

If  we  admit  weak  solutions  into  the  discussion  of  a problem, 
(such  as  those  of  the  type  introduced  in  Sections  4.0  - 6.0),  we  would 
anticipate  that  since  the  admissible  class  of  solutions  has  been  greatly 
widened,  there  could  possibly  be  many  solutions  to  that  problem.  It 
is  well  known  that  this  is  indeed  the  case  in  the  theory  of  quasi-linear 
hyperbolic  partial  differential  equations.  See  for  example  Lax  [8]. 
The  boundary-value  problem  referred  to  at  the  end  of  the  preceding 
section  confirms  this  to  be  the  case  in  the  present  context  as  well. 

In  such  circumstances,  it  is  essential  to  introduce  criteria 
which  single  out  a physically  admissible  solution  from  among  the  many 
solutions  admitted  by  the  differential  equations.  The  second  law  of 
thermodynamics  appears  to  play  such  a role  in  gas  dynamics.  Lax  [8] 
has  examined  "entropy  conditions"  which  furnish  such  criteria  in  the 
context  of  hyperbolic  systems  of  conservation  laws. 

An  analogous  "entropy  condition"  in  the  context  of  elastostatics 
was  proposed  by  Knowles  and  Sternberg  [3]  and  subsequently  extended 
by  Knowles  [4],  A thermodynamic  motivation  for  the  proposed  condi- 
tion, in  the  case  of  compressible  materials,  was  also  given  in  [4].  In 
the  three-dimensional  case,  a quasi-static  time  dependent  family  of 
equilibrium  states  was  considered,  the  time  merely  playing  the  role 
of  a history  parameter,  and  it  was  then  required  that 


L 
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JsN-vdA--^  Jw(F)dv^O  (7.  1) 1 

d£  & 

for  every  regular  sub-domain  £ of  ft  , at  each  instant  of  the  time  in- 
terval considered.  Here  t is  the  time  and  v the  quasi- static  particle 
velocity.  Equation  (7.1)  gives  expression  to  the  idea  that  the  rate  at 
which  elastic  energy  is  being  stored  in  & cannot  exceed  the  rate  at 
which  work  is  being  done  on  £ . 

One  shows  easily  that  for  a sub-domain  £ of  the  body  which  is 
such  .hat  the  field  quantities  have  classical  smoothness  properties  at 
each  interior  point,  the  global  condition  (7.  1)  holds  with  inequality 
replaced  by  equality  by  virtue  of  the  field  equations.  This  is  indeed  as 
one  would  expect,  and  accordingly  (7.  1)  imposes  no  local  restrictions 
at  a point  where  the  fields  are  smooth.  If  however  an  elastostatic 
shock  is  present  in  the  domain  fi  , then  (7.  1)  does  not  hold  automati- 
cally and  consequently,  at  each  point  on  the  shock  it  imposes  a local 
restriction  on  the  jumps  of  the  field  quantities. 

Now  consider  a quasi- static  family  of  plane  strain  piecewise 
homogeneous  elastostatic  shocks  in  a homogeneous,  incompressible, 
elastic  solid.  It  can  be  shown  that,  if  at  some  instant  t (7.  1)  holds 
with  strict  inequality  for  all  sub-domains  £ which  intersect  at  , then 

(i)  the  motion  of  the  shock-line  at  that  instant  is  translatory 
in  a direction  not  parallel  to  itself.  Moreover,  if  we  orient 

the  shock-line  at  such  that  this  translation  is  directed  into 

+ 

II  , then  at  that  instant 

*Body  forces  were  omitted  from  this  discussion. 
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tW(F)-OopNpFovNv]>0  . (7.2) 

Conversely,  if  at  some  time  t the  quasi-static  family  of  solutions  con- 
forms with  (i),  then  (7.  1)  holds  with  strict  inequality  at  that  instant. 

On  the  other  hand  we  can  show  that  if  at  the  instant  t (7.  1) 
holds  with  equality  for  all  sub-domains  & then  either 

(ii)  the  shock-line  is  instantaneously  stationary  at  that 
moment, 

or  (iii)  the  shock-line  s£  is  instantaneously  in  a state  of  translation 
parallel  to  itself  at  that  moment, 

A 

or  (iv)  the  jump  of  W(F)  - o^N^F^N^  across  the  shock  is  zero, 

(in  which  case  the  shock-line  motion  is  not  restricted  to 
being  translatory). 

Conversely,  if  at  some  time  t the  quasi-static  family  of  solutions  con- 
forms with  one  of  (ii),  (iii),  and  (iv),  then  (7.  1)  holds  with  equality  at 
that  instant. 

Finally  one  can  show  that  if  at  some  time  t (7.  1)  holds  for  all  sub- 
domains  & , and  if  in  addition  it  holds  with  equality  for  some  sub-domain 
which  intersects  the  shock,  then  in  fact,  at  that  instant  (7.  1)  holds  with 
equality  for  all  sub-domains  & . We  conclude  from  this  that  the  preceding 
are  the  only  possibilities.  Therefore,  if  (7.  1)  holds  it  is  necessary 
that  one  of  (i)  - (iv)  hold.  Conversely  if  one  of  (i)  - (iv)  holds  this  is 
sufficient  to  ensure  that  (7.  1)  hold. 

One  arrives  at  (i)  - (iv)  by  applying  to  the  incompressible  case 
the  parallel  arguments  used  by  Knowles  and  Sternberg  in  [3],  or  by 
specializing  to  this  context  the  results  of  Knowles  [4].  Since  (7.  1) 
implies  that  the  presence  of  an  elastostatic  shock  decreases,  or  at 
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least  does  not  increase,  the  stored  energy  in  the  body,  we  refer  to 
(7.  2)  as  the  dissipativity  inequality. 

It  is  apparent  from  (i)  - (iv)  that  the  dissipativity  requirement 
(7.  1)  may  be  viewed  as  restricting  the  admissible  class  of  quasi-static 
motions.  The  only  quasi-static  motions  admitted  by  it  are  those  in 

the  value  of  { - W (F)  -t-  FN*ON}  at  a particle  does  not  decrease  as 
the  particle  crosses  the  shock-line. 

It  may  be  remarked  that  the  dissipativity  inequality  does  not 
rule  out  any  piecewise  homogeneous  elastostatic  shock  itself  as  being 
inadmissible,  since  any  given  elastostatic  shock  can  always  be  embedded 
*n  a suitable  time  dependent  family  of  such  shocks  which  conforms  with 
the  dissipativity  inequality. 

As  one  would  expect,  and  as  is  verified  by  Knowles  [4],  these 
results  remain  true  locally  in  the  general  case  of  a curved  shock  in  a 
non-homogeneous  elastic  field,  with  the  exception  that  the  shock  motion 
may  no  longer  be  restricted  to  translation.  The  latter  property  is 
clearly  peculiar  to  piecewise  homogeneous  elastostatic  shocks. 

Using  (2.4),  (2.6),  (4.  13),  (5.3),  and  evaluating  the  left  hand 
side  of  (7.2)  in  the  frame  X'  leads  to 

fW(F)-FftpOaYNpNY]+=[W(F)]+  + KT^  , (7.3) 

where  we  have  also  used  the  fact  that  Tj£  is  continuous  across  • 
Therefore  the  inequality  (7.2)  may  be  written  in  the  simpler  form 

[W(F)]++kt*'>0  . (7.4) 

where  £ is  presumed  to  be  oriented  such  that  it  moves  into  il-  as  time 
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t increases. 

In  the  particular  case  when  the  material  at  hand  is  isotropic, 
we  have  from  (2.33),  (6.  10)  and  (6.  18)  that 


y ' + 

T?.  = 2bW'(I)  . 


(7.5) 


whence  (7.4),  by  virtue  of  (6.  14),  (6.  18)  and  (6.  19),  may  be  equivalently 


written  as 


W(I)  - W(I  + 2bK  +CK2)  + 2bHW,(I)>0  . 


(7.6) 


Note,  however,  that  (6.23)  may  alternatively  be  written  as 


1 flW  + 7 + 

Mx)  = (I  + 2bx  + cx^)  - bW#(I ) , 


(7.7) 


whence  (7.6)  takes  the  simple  form 


J h(x)dx<  0 


(7.8) 


We  will  make  use  of  this  form  of  the  dissipativity  inequality  in  the 
example  taken  up  in  the  next  section. 

Finally,  we  return  to  anisotropic,  incompressible,  elastic 

solids  in  order  to  determine  the  weak  shock  approximation  to  the  value 

of  the  jump  of  fW(F)-Fa(3OavNpNv]  across  the  shock -line.  Recall 

from  Section  5.  1,  where  we  first  looked  at  weak  elastostatic  shocks, 

that  we  now  assume  that,  given  the  deformation  gradient  F with  unit 

determinant  and  the  pressure  p , there  exist  functions  0(h)  and  p(n), 

both  sufficiently  smooth  in  a neighborhood  of  * =0,  such  that  F (k) 

orp 


1 
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defined  by 


Fop(K)=^crp+Klcr(H)nY(H)^Yp  * (7*9) 

conforms  with  the  traction  continuity  requirement  (4.  12).  Observe 
from  (7.  9)  that 


^P<K>=*:Jp-K,Y<KV,)*:«'v  ' <7-‘01 

It  is  first  necessary  to  analyze  the  traction  continuity  condition 
(4.  12).  To  this  end  set 

^(H)=JopNp(K)-aap(K)Np(K)  . (7.11) 

which  because  of  (4.9),  (5.3),  (7.  10)  and  the  perpendicularity  of  ^ 
and  n leads  to 


r8W(F)  8W(F(k)) 


8F 


8F_ 


op  op 


(7.12) 


Differentiating  (7.  12)  with  respect  to  k and  using  (2.  11),  (5.  10)  and 
(7.9)  gives 


= " caPY6^^ir6^p^0^Y^0^nfl^0^  + **p(0)  ' <7'13» 

The  continuity  of  traction  across  the  shock  requires  that 

Aa(K)  = 0 (7.14) 

for  all  sufficiently  small  h , from  which  it  follows  that  in  particular 

A;(0)  = 0 . 


(7. 15) 
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From  (7.  13)  and  (7.  15)  we  find  that 

C«Pv6(E)^6'v(0»n,(0)Np(0|-P,(0^|ivlNp(0,  = 0 • ,7'16’ 

As  one  would  expect,  (7.  16)  is  in  fact  the  same  as  (5.  12)  because  of 

(5.3) .  Differentiation  of  (7.  12)  twice  with  respect  to  * , together  with 
the  symmetry  capy6=  cY6c>p  » the  fact  that  ~*H  = 0 and  (2**1)*  (5.1)  - 

(5.3) ,  (5.10)  and  (7.9)  leads  to 

V0)Aa<°>  = ' c2<°^Pv»M(E,l«(0)iV(0)iX  (0)Np(0)N6(0)Nu(0) 

-4c<0>capY6(£)<a(0)lY(0)V°)N6(0) 

-2p'(0)FpY1Np(0)f/fiy(0)  + 2p'(0)F(^1N/p(0)fa(0)  , (7.17) 

where  we  have  set 

83W(F) 

WM<P-»r^ry,lr>M  - (7‘ 181 

and  c(h)  was  defined  in  (5.  1).  Because  of  (7.  14)  we  have  that 
A*(0)  =0  , whence  we  have  from  (7.  17)  that 

4c<0)c«PY6(^i«(°)iV(0)NP(0)Ni(0) 

= -c2<0)d«PY6)iu^^(0)<Y<0)i)i(0)Np(0)N6(0)Nu(0) 

-2p/(0)FpY1Np(0)/;(0)  + 2p'(0)Fpcr1N'p(0)fa(0)  . (7.19) 

A 

We  now  compute  the  jump  in  { W(F)  - FN*£N}  across  the  shock. 
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To  this  end,  let 

n<» ) = W(F(« ))  - W(F)  + FapNp(H  )JavNy(K ) 

-Fap(-'Np(H)^v(-)NY(«)  . (7.20) 

Because  of  traction  continuity,  the  fact  that  ,i*n  = 0 , (4.9),  (5.3)  and 
(7.9)  we  can  write  (7.20)  as 

/v  «• 

ti(h)  = W(F(h))- W(F)-HC(H)fa(H)N  (H)|^^l)  . (7.21) 

Y r ay 

Clearly, 

11(0)  =0  , (7.22) 

by  virtue  of  (7.9).  Differentiating  (7.  21)  with  respect  to  k and  using 
(2.11),  (5.  3)  and  (7.  9)  gives 

Tl/(H)  = -Hc(H)C(>pY6(F(K))iY(K)N6(K)|r{Kc(H)i(y(H)N(3(H)}  , (7.23) 

from  which  we  have  that 

T)'(0)  = 0 . (7.24) 

Differentiating  (7.  23)  with  respect  to  k and  using  (7.9)  leads  to 

T1*(0)  = - c2(0)CofpY6(F)fa(0)f  Y(0)Np(0)N6(0)  , (7.  25) 

which  because  of  (5.3)  and  (7.16)  gives 


I 


Ti*(0)  = -p'(0)na(0)<a(0)  , 


(7.26) 
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which  in  turn,  because  l Q(0)na(0)  = 0 , implies  that 


h"(0)  = 0 


(7.27) 


Finally,  differentiating  (7.23)  twice  with  respect  to  h , using  the  sym- 
metrV  copY6=Cv6a(3  ’ (5*3)'  (7.9)  and  (7.  18)  leads  to 


T]m(0)  = - 2c  (0)d 


aPySXM^V^V0^  (0)Np(0)N6(0)NM(0) 


’ 6c  (0)caPY6(E)io(0)\(0)N6(0)Np(°) 

'6c(0)copY6(^MY(0)Np(0)N6(0)fc/(0)lo(0)  + c(0)ii(0)}  (7.28) 

which  on  using  (5.3),  (7.  19)  and  f ^(OJn^O)  = 0 implies  that 


rT(0)  = - ^c  (0)d 


2C  (0>d<rpY6Xu(^,o(0)1Y(0)iX(0>V0)N6{0)NM(0) 


-3p'(0)[n'a(0)fo(0)+no(0)<;(0)} 


Since  the  vectors  ^(k)  and  n(n)  are  perpendicular  to  each  other. 


(7.29) 


fa(*Ono(K)  = 0 for  all  sufficiently  small  k 


(7.30) 


Differentiating  (7.30)  with  respect  to  h shows  that 


i;(0)no(0)  + y0)n^(0)  = 0 , 


(7.31) 


so  that  finally  we  may  write  (7.29)  as 


n"(0>  = -2C  <°)dapv6XM(5:)V0)iY(0)iX(0)V0)N6(0)NU(0)  * (7'32) 
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Therefore  (7.21),  (7.22),  (7.  24)  and  (7.  27)  allow  us  to  write 

^(£)-FapVavNvJ+  = iT1"(0)K3+o(K3)  as  *-°  • <7-33> 

where  T)w(0)  is  given  by  (7.32).  We  observe  that  the  jump  in 

fW(F)-  FoP^(3aoyNy^  across  shock  is  of  the  third  order  in  the  shock 

strength  x , which  is  as  in  the  case  of  compressible  elastic  solids. 

This  is  analogous  to  the  situation  in  gas  dynamics  where  the  entropy 
jump  is  of  the  third  order  in  the  appropriate  shock  strength.  * 


See  references  cited  in 


[4]. 
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f 


8.  1 An  Illustrative  Example 

For  the  purpose  of  illustrating  the  results  of  the  previous  sec- 
tions and  demonstrating  how  in  a particular  case  one  could  in  fact  ob- 
tain even  more  information  than  has  been  indicated,  we  now  specialize 
our  constitutive  law.  Consider  the  hypothetical  class  of  homogeneous, 
isotropic,  incompres sible,  elastic  solids  for  which  the  plane  strain 
elastic  potential  is  given  by 

W(I)  = |akg{l  -exp(-  , H>0  , kQ>0  . (8.1) 

2ko 

One  sees  immediately  from  (8.  1)  that  (2.40)  is  satisfied  whence  this 
class  of  materials  has  a positive  shear  modulus. 

According  to  (3.32),  we  have  in  simple  shear,  the  shear  stress- 
amount  of  shear  relation 


T(k)  =|ikexp^ — ^“2)  • 


(8.2) 


A sketch  of  the  response  curve  in  shear  defined  by  (8.2)  is  shown  in 
Fig.  2.  The  significant  feature  of  this  for  our  purposes  is  that  T'(k) 
is  positive  for  all  k in  the  interval  (-  kg,  kg)  and  is  non-positive 
otherwise.  The  implications  of  this  as  far  as  the  issue  of  the  ellipti- 
city  of  the  displacement  equations  of  equilibrium  are  concerned  were 
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observed  in  Section  3. 

We  now  turn  to  the  issue  of  piecewise  homogeneous  elastostatic 

shocks.  Suppose  we  are  given  the  deformation  gradient  F , and  hence 

(through  (2.4),  (6.5)  and  (6.9))  the  associated  value  of  p (say  (3)  , and 

+ + 

the  pressure  p on  n . We  look  at  the  question  of  the  existence  of  a 

A 

corresponding  elastostatic  shock  with  spatial  shock-line  inclination  0 

A A 

to  the  yj-axis.  0 and  p are  held  fixed  in  this  discussion,  and  as 

A A 0 

noted  previously  we  may  assume  (0  , P)  to  be  in  G , with  no  loss  of 
generality.  We  recall  that  a corresponding  piecewise  homogeneous 
elastostatic  shock  exists  if  and  only  if  the  function 


h(x)  = (b  + c x)W'(?+  2bx  + cx^)  - bW'(?)  , 


(8.3) 


where 


A A 


J l-P* 


sin  20  <0  , c = 1 -PPos-±&  >o  , 


J 1-P2 


+ ? 

I = ---  >2  , 


(8.4) 


V 1-]^ 


has  a zero  at  some  x^O.  Using  (8.1)  in  (8.3),  we  find,  for  the  type 
of  materials  under  consideration,  that 


w * ^){(b*cx).xp(- i2n&±~h) . b}  . 


(8.5) 


Case  (i)  Suppose  0 and  p are  such  that 


See  discussion  following  Equation  (3.33). 


-68- 


b2>cko  * (8.6) 

Then  h/(0)<0  . One  shows  easily  from  (8.4),  (8.5)  and  (8.6)  that  in 

this  case  h(x)  has  a unique  zero  (in  addition  to  the  one  at  the  origin) 

at  x = k , where  k is  a positive  number  and  is  such  that  Jh(x)dx<0  . 

0 

It  follows  that,  corresponding  to  the  homogeneous  deformation 

A A 

associated  with  (3  on  II  and  to  the  inclination  0 compatible  with 

(8.4)  and  (8.6),  there  exists  a unique  piecewise  homogeneous  elasto- 
static  shock  with  positive  shock  strength  k . Furthermore,  suppose 
this  piecewise  homogeneous  shock  is  embedded  in  a quasi- static  family 
of  shocks.  Then  if  at  the  instant  when  the  family  of  shocks  coincides 
with  this  given  shock  the  shock-line  £ is  translating  into  II  , it  con- 
forms with  the  dissipativity  inequality. 

A A 

Case  (ii)  Suppose  0 and  p are  such  that 

b2<ckQ  . (8.7) 

Then  h/(0)>0  . In  this  case,  it  is  easily  verified  by  virtue  of  (8.4), 

(8.5)  and  (8.7)  that  h(x)  has  a unique  zero  (in  addition  to  the  one  at 
the  origin)  at  x = k , where  k is  a negative  number  such  that 

K 

J h(x)dx>0  . 

0 

It  follows  that  corresponding  to  the  homogeneous  deformation 

A A 

associated  with  P on  n and  the  inclination  0 compatible  with  (8. 4) 
and  (8.7),  there  exists  a unique  piecewise  homogeneous  elastostatic 
shock  with  negative  shock  strength  h . Furthermore,  suppose  this 
piecewise  homogeneous  shock  is  embedded  in  a quasi- static  family  of 
shocks.  Then,  if  at  the  instant  when  the  family  of  shocks  coincides 
with  this  given  shock  the  shock-line  ^ is  translating  into  II  , it 
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conforms  with  the  dissipativity  inequality. 

A A 

Case  (iii)  Suppose  0 and  p are  such  that 


b 


2 


(8.8) 


Then  h'(0)  = 0.  In  this  case  the  only  zero  of  h(x)  is  at  the  origin, 

*4* 

from  which  we  conclude  that  if  the  homogeneous  deformation  on  II  is 

A /\ 

such  that  the  associated  value  of  P and  the  (proposed)  inclination  0 
conform  with  (8,4)  and  (8.8),  then  there  is  no  corresponding  piecewise 
homogeneous  elastostatic  shock. 

These  results  are  best  visualized  on  the  (0,p)-plane.  Using 
(8.4)  we  have  that 


2 - p2cos2  20  + k2p >/ 1 ~ p2  cos20+(p2-ko2^T7) 

ck°=  ITFI  ' 


(8.9) 


Let  T be  the  curve  in  the  first  quandrant  of  the  0 - P plane  whose 
equation  is 


T : p2cos220  - k2p 7 1 - p2  cos  20  - (p2-  k ZQJ 1 - p*)  = 0 . (8.10) 


o 

T separates  G into  two  regions  as  shown  in  Fig.  3.  Case  (i) 
refers  to  points  in  the  hatched  open  region  jhown  there,  while  Case  (ii) 
refers  to  points  in  the  unhatched  open  region.  Points  on  T refer  to 
Case  (iii).  One  finds  that  T has  a minimum  point  at  (0  , P ) where 

C 0 


0 - -r  CO  S 
e 2 


ko^ko*4 

k2  + 2 


(8. 11) 
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From  (2  43 'i  (".33),  (8.2)  and  (8.4)  we  find  that  the  displacement 

+ 

equations  of  equilibrium  are  elliptic  on  II  , if  and  only  if  the  deforma- 
tion there  is  such  that  the  associated  value  of  p is  less  than  P . 

+ 

Suppose  that  the  given  deformation  on  II  is  such  that  the  displacement 

equations  of  equilibrium  are  non-elliptic  there.  Then  • The 

spatial  characteristics  associated  with  this  deformation  are  inclined 

+ 

to  the  Xj-principal  axis  of  G at  angles  or  , which  because  of  (3.46), 
(8.1)  and  (8.4)  are  given  by 


cos2o=  (1 

2P 

Note  however,  that  the  equation  of  the  curve  T , (8.  10),  can  alterna- 
tively be  written  as 


f : cos  20 


k^l-p2±({(kQ  + 2)Vl-P^-2j  {(k2-2)Vl-p2-2})i 
= — 


(8.13) 


It  is  immediately  evident  from  a comparison  of  (8.12)  and  (8.  13)  that, 
the  abscissa  of  the  points  on  T corresponding  to  P^P  give  the  spa- 
tial  characteristic  inclinations  corresponding  to  the  deformation  asso- 

A 

ciated  with  p . 

We  now  summarize  our  findings  for  the  particular  class  of 

materials  at  hand.  Corresponding  to  any  given  homogeneous  deforma- 

+ 

tion  on  II  we  can  have  a piecewise  homogeneous  elastostatic  shock 

+ f.  , 

(provided  F is  not  proper  orthogonal,  i.e.  PrO). 

If,  at  the  given  deformation,  the  displacement  equations  of 

+ « 

equilibrium  are  elliptic  on  II  , so  that  P<Pe  » the  spatial  shock-line 
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may  be  inclined  at  any  angle  0 provided  it  is  not  parallel  to  the  princi- 
pal axes  of  G (i.e.  0^0  , . One  can  show  that  for  such  an  elasto- 

static  shock,  the  displacement  equations  of  equilibrium  are  non-elliptic 
on  II  . Furthermore  the  corresponding  shock  strength  is  negative  and 
a quasi-static  motion  from  such  a configuration  is  compatible  with  the 
dissipativity  inequality  if  the  shock  moves  into  II. 

On  the  other  hand  if  the  displacement  equations  of  equilibrium 

+ A 

are  non-elliptic  at  the  given  deformation  on  II  , so  that  p^p^  , the 

spatial  shock-line  may  be  inclined  at  any  angle  0 provided  it  is  not 

+ 

parallel  to  the  principal  axes  of  G nor  parallel  to  the  spatial  charac- 
teristic directions  associated  with  the  deformation  on  S . In  this  case 
the  sign  of  the  shock  strength  and  the  admissible  direction  of  quasi- 
static motion  depends  on  the  specific  shock-line  inclination  (see  Fig.  3). 
In  particular  note  that  the  admissible  direction  of  quasi-static  shock 
motion,  for  dissipativity,  is  governed  solely  by  whether  the  spatial 
shock-line  inclination  is  between  or  outside  the  inclinations  of  the 

2 spatial  characteristics  (in  the  relevant  quadrant)  associated  with  the 

+ 

deformation  on  n . The  ellipticity  or  non-ellipticity  of  the  displace- 
ment equations  of  equilibrium  on  n also  turn  out  to  depend  on  the 
specific  shock-line  inclination. 

Acknowledgment 

The  author  wishes  to  express  his  sincere  appreciation  to 
Professor  James  K.  Knowles  for  his  invaluable  guidance  throughout 
this  investigation. 


-72- 


REFERENCES 


[1]  J.  K.  Knowles  and  Eli  Sternberg,  On  the  ellipticity  of  the  equations 
of  nonlinear  elastostatics  for  a special  material.  Journal  of 
Elasticity,  5 (1975),  3-4,  p.341. 

[2]  J.K.  Knowles  and  Eli  Sternberg,  On  the  failure  of  ellipticity  of  the 
equations  for  finite  elastostatic  plane  strain.  Archive  for  Rational 
Mechanics  and  Analysis,  £3  (1977),  4,  p.321. 

[3]  J.K.  Knowles  and  Eli  Sternberg,  On  the  failure  of  ellipticity  and  the 
emergence  of  discontinuous  deformation  gradients  in  plane  finite 
elastostatics.  Journal  oi  Elasticity,  8 (1978),  4,  p.329. 

[4]  J.K.  Knowles,  On  the  dissipation  associated  wifh  equilibrium  shocks 
in  finite  elasticity.  Journal  of  Elasticity,  9 (1979),  2. 

[5]  J.R.Rice,  The  localization  of  plastic  deformation.  Proceedings  of 
the  14th IUT AM  Congress,  (W.T.Koiter,  editor).  Delft,  Netherlands, 
August  1976,  p.207. 

[6]  C.Truesdell  and  W.Noll,  The  nonlinear  field  theories  of  mechanics, 

. Handbuch  der  Physik,  Vol.III/1,  Springer,  Berlin,  1965. 

[7]  P.  Chadwick,  Continuum  Mechanics,  Wiley,  New  York,  1976. 

[8]  P.D.Lax,  Hyperbolic  systems  of  conservation  laws  and  the  math- 
ematical theory  of  shock  waves.  Regional  Conference  Series  in 
Mathematics,  No.  11,  Society  for  Industrial  and  Applied  Mathe- 
matics, Philadelphia,  1973. 


ADMISSIBLE  SHOCK  MOTION  IS  INTO  II 


| | ADMISSIBLE  SHOCK  MOTION  IS 


FIGURE  3.  PLANE  OF  PARAMETERS 
ADMISSIBLE  REGION  CL 


INTO  n 


4>  AND  £ 


SECURITY  CLASSIFICATION  OF  This  PAGE  (When  Dale  Entered] 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1 REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

No. 41  ^ 

3.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  (end  Subtitle) 

Discontinuous  deformation  gradients  in  plane 
finite  elastostatics  of  incompressible  materials 

5.  TYPE  Of  REPORT  A PERIOD  COVERED 

Research 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORfsJ 

Rohan  C.  Abeyaratne 

8.  CONTRACT  OR  GRANT  NUMBERO) 

N00014-75-C-0196  ✓ 

».  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

California  Institute  of  Technology 

Pasadena,  California  91125 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  6 WORK  UNIT  NUMBERS 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 

Washington,  D.  C. 

12.  REPORT  DATE 

May,  1979 

13.  NUMBER  OF  PAGES 

76 

14.  MONITORING  AGENCY  NAME  S ADDRESSfl t dl Iterant  Irom  Controlling  Oltlee) 

15.  SECURITY  CLASS,  (ol  thle  report) 

Unclassified 

1S«.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (ot  thle  Raport) 

Approved  for  public  release;  distribution  unlimited 

17.  DISTRIBUTION  STATEMENT  (ol  the  abatrael  entered  In  Block  30,  II  dlllerent  trow  Report) 

16.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Continue  on  revaraa  aide  II  naceeaary  and  Identity  by  block  number) 

Elastostatics 

Incompres  sible 
! Nonlinear 

Discontinuous  deformation  gradients 

2oV  ABSTRACT  (Continue  on  reveree  tide  If  neceeemry  end  Identity  by  block  number)  J 

uThis  investigation  is  concerned  with  the  possibility  of  the  change  of  type  of 
fVio  riifft* -rpnti al  equations  governing  finite  plane  elastostatics  for  incompres- 
sible  elastic  materials,  and  the  related  issue  of  the  existence  of  equilibrium 
fields  witn  discontinuous  deformation  gradients.  Explicit  necessary  and 
sufficient  conditions  on  the  deformation  invariants  and  the  material  for  the 
ellipticity  of  the  plane  displacement  equations  of  equilibrium  are  establishec  . 

DD  ,^nRM73  1473  EDITION  OF  I NOV  6S  IS  OBSOLETE 
S/N  0102  LF  014  6601 


SECURITY  CLASSIFICATION  OF  TNI! 


(When  Data  Entered) 


$ 


unclassified 


SECURITY  CLASSIFICATION  OF  this  PAGEfNftOT  Dm tm  Kntmrmil) 


The  issue  of  the  existence,  locally,  of  "elastostatic  shocks"  — elastostatic 
fields  with  continuous  displacements  and  discontinuous  deformation  gradient 
is  then  investigated.  It  is  shown  that  an  elastostatic  shock  exists  only  if  the 
governing  field  equations  suffer  a loss  of  ellipticity  at  some  deformation. 
Conversely,  if  the  governing  field  equations  have  lost  ellipticity  at  a given 
deformation  at  some  point,  an  elastostatic  shock  can  exist,  locally,  at  that 
point.  The  results  obtained  are  valid  for  an  arbitrary  homogeneous,  iso- 
tropic, incompressible,  elastic  material. 


JiruSITV  CLASSIFICATION  OF  THIS  ^AOtfWiFo  Knfmtfl 


